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^ , Abstract 

bX)' In this (expository) paper and its (more technical) companion we reformu- 

^ , late some results of the Nobel Prize winning paper by Werner Heisenberg into 

' modem mathematical language of honeycombs. This language was recently de- 

, veloped in connection with complete solution of the Horn problem (to be defined 

CSJ ' and explained in the text). Such a reformulation is done with the purpose of 

posing and solving the following problem: is by analyzing the (spectroscopic) 
experimental data it possible to recreate the underlying microscopie model gen- 
erating these data? Although in the case of Hydrogen atom positive answer is 
known, to obtain an afhrmative answer for spectra of other quantum mechan- 
ical Systems is much harder task. Development of Heisenberg's ideas happens 
OO ' to be the most useful for this purpose. It supplies needed tools to solve the 

Veneziano and other puzzles. The meaning of the word "puzzle" is two-fold. 
On one band, it means (in the case of Veneziano amplitudes) to find a physical 
model reproducing these amplitudes. On another, from the point of view of 
combinatorics of honeycombs, it means to find explicitly fusion rulcs for such 
amplitudes. Solution of these tasks is facilitated by our earlier developed string- 
theoretic formalism. In this paper only qualitative arguments are presented 
, (with few exceptions). These arguments provide enough evidence that the un- 

derlying model compatible with Veneziano amplitudes is the standard (i.e. non 
supersymmetric!) QCD. In addition, usefulness of the proposed formalism is 
■ illustrated on numerous examples such as physically motivated solution of the 

saturation conjecture (to be defined in the text), derivation of the Yang-Baxter 
and Knizhnik-Zamolodchikov equations,Verlinde and Hecke algebras, computa- 
tion of the Gromov-Witten invariants for small quantum cohomology ring, etc. 
Finally, we discuss possible uses of these ideas in condensed matter physics 
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1 Introduction 



1.1 Some general facts about Veneziano condition and its 
CFT analog 

To avoid repetitions, we refer our readers to earlier papers, Refs.[l-3], which 
shall be called Parts I-III respectiveljQ. In particular, in Part I we noticed tliat 
for the 4-particle scattering amplitude the Veneziano condition is given by 

a{s) + a{t) + a{u) ^ -1, (1.1) 

where a(s), a(t),a{u) e Z . Eq.(l.l) is a simple statement about the energy - 
momentum conservation. Although the numerical entries in this equation can be 
changed to make them more suitable for theoretical treatments, actual physical 
values can be subsequently re obtained by appropriate coordinate shift. Such 
a procedure is not applicable to amplitudes in conformai field theories (CFT) 
where the periodic (antiperiodic, etc.) boundary conditions cause energy and 
momenta to become a quasi -energy and a quasi- momenta in terminology taken 
from the solid state physics. Noticed differences between the combinatorial 
properties of CFT and high energy physics amplitudes the major theme of this 
work. 

To explain things better, we would like to rewrite Eq.(l.l) in a more con- 
venient form. FoUowing Ref. [4], without loss of generali ty, the honiogenous 
equation 

a{s)m + a{t)n + a{u)l + k ■ l = Q (1.2) 

where m, n, l, k are some integers can be added to Eq.(l.l) thus producing the 
Veneziano-type equation: 

a(s)m + a{t)n + a{u)l — k. (Veneziano) (1.3a) 

This equation is equivalent to 

ni+n2+n3 = iV, (1.3b) 

where hy design ali entries are nonnegative integers. For the case of multiparticle 
scattering we anticipate that this equation is to be replaced by 

no + --- + nfc=iV (1.4) 

as discussed in Part II. Combinatorially, the task lies in finding ali nonnega- 
tive integer combinations of no,...,nfc satisfying Eq.(1.4). It should be noted 
that such a task makes sense as long as N is assigned. But the actual value 
of N is not fixed and, hence, can be chosen quite arbitrarily. As explained 
in Part I, the value of N should coincide with the exponent of the Fermat 
(hyper)surface if (in contrast with traditional string-theoretic treatments) we 

^In this work, when we occasionally refer to equations in other parts, we shall use the 
following convention, e.g. Eq.(I 3.28), refers to Eq.(3.28) of Part I, etc. 
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interpret Veneziano amplitudes as periods of the Fermat (hyper) surfaces living 
in the complex projective space. This requirement is not too rigid, however, 
in view of Eq. (1.3.29). Indeed, the mathematica! statement of the type given 
by Eq.(1.4) should be considered before the bracket operation <...> defined in 
Part I is appUed. This means that we shaU be working mainly with the pre- 
cursors of the period integrals and, therefore, we can choose any non negative 
integer value for A^.Physically correct value of N can then be reobtained by the 
appropriate coordinate shift. 

In CFT such shift looses its meaning due to periodicity. In this case since the 
Veneziano condition, Eq.(1.3a), should be replaced by the Kac-Moody-Bloch- 
Bragg (K-M-B-B) condition (e.g. see Eq. (1.3.22)): 

a{s)m + a{t)fi + a{u)ì^mN + nN + lN, {K ~ M ~ B - B) (1.5) 

where N is the same as before but m, n, l are arbitrary integerfl This circum- 
stance causes the energy (in free space) to become a quasi-energy (in sohds). 
Eq.(1.5) is the most fundamental equation in X-ray crystallography [5] where 
it is known as the Bragg condition. In sohd state physics essentially the same 
equation is known as the Bloch equation. As it foUows from monographs by Kac 
[7] and, much earlier, by Bourbaki [8], the affine Lio groups and the associated 
with them Lie algebras are generalizations of the Weyl-Coxeter reflection groups 
made by analogy with crystallographic groups in solid state physics. In the lan- 
guage of solid state physics, the Weyl-Coxeter reflection groups are "point" 
groups while the affine groups are made of semidirect products of translation 
groups and point groups are "spatial" groups. Unlike the solid state physics, in 
the present case translations can be performed in the Euclidean, hyperbolic or 
spherical spaces. These spaces need not be 3 dimensionai. To make a comparison 
with the existing literature, we shall be concerned only with the Euclidean-type 
translations. For a quick concise review of ali these concepts we refer our readers 
to the Appendix A of Part II. 

The arbitrariness of choosing N in Eq.(1.4) represents a kind of gauge free- 
dom. As in gauge theories, we can fix the gauge by using some physical con- 
siderations. These include, for example, an observation made in Part I that 
the 4-particle amplitude is zero if any two entries info Eq.(l.l) (or, which is 
the same, into Eq.(1.3b)) are the same. This fact prompts us to arrange the 
entries in Eq.(1.3b) in accordance with their magnitudes, e.g. ni > ri2 > ris. 
More generally (in view of Eq.(1.4)), we can write: no > ni > • • • > rifc > 
Provided that Eq.(1.4) holds, we shall cali such a sequence a partition and de- 
note it as A=(no, .■■,nk)- If A is a partition of A^,then we shall write X\- N. It 
is well known from combinatorics [4,9] that there is one-to-one correspondence 
between the Young diagrams and partitions. We used this observation in Part II 
for designing new partition function capable of rcproducing the Veneziano (and 

^Clearly we can combine them together but we do not do this to keep an analogy with the 
solid state physics. 

•^The last inequality: nj. > 1, is chosen only for the sake of comparison with the existing 
literature conventions, e.g. see Ref. [6]. 
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Veneziano-like) amplitudes. Now we would like to look at the same problem 
from a difFerent angle. 

1.2 The additive Horn problem 

Consider some k x k Hermitian matrix Ti. whose spectrum A = {Ai > ... > Afe) 
can be written as a partition made of weakly decreasing sequence of k reao 
numbers. Conversely, for every spectrum A there is a set Ox of Hermitian 
matrices whose spectrum is A. Suppose that we are given 3 such spectra: A, 
and z/, then there should be matrices Hx G 0\, Ti^ G and H,^ E Ov such that 
T-l\ + = Tiy. More precisely, according to Hermann Weyl [10], the foUowing 
problem can be formulated. 

Problem 1.1. Assuming the eigenvalues of two k x k Hermitian matrices 
T-Lx and Ti^ to be known, how does one determine ali possible sets of eigenvalues 
for the sum Hx + U^M 

For k = 1 the answer is obvious but for k > 1 the answer is much less 
obvious. Since tr[T-lx + Ti^] — trpi.^] we always bave the trace condition: 

ui + ... + Uk = Al + ... + Afe + /il + ... + ^Jik■ (1-6) 

Clearly, in addition, we can expect that vi < A1+/I1 since the largest eigenvalue 
of TCx + Ti.^ is at most the sum of TCx and individuai eigenvalues. 

Let now I, J and K be subsets of {1, k} then, what can be said about the 
validity of the inequality 

$^^.<EA, + ^M. ? (1-7) 
iei jeJ fceK 

Such type of inequalities were analyzed by Horn [13] who formulated the fol- 
lowing. 

Conjecture 1.2. (Horn) A triple A, /i and ly can represent the eigenvalues of 
the k X k Hermitian matrices Hx , and Tii, , where Ti.x + Ti. fi = Ti-v , if and only 
if the trace condition, Eq.(1.6), holds and, in addition, the inequality given by 
Eq.(1.7) holds for ali sets {I, J and K} G , r < k. is defined recursively 
as foUows. For each positive integer k and r < k, let 

U::^{{I,J,K)\Y,^ + Y,J=J2k + r{r + ^)/2}. (1.8) 
iei jeJ fceK 

^In accord with the existing mathematica! hterature we refer to this (the most general 
case) as "classical". Accordingly, the "quantum" case corresponds to a situation when ali 
numbers are integers. We shall say more on this topic in Section 5. 

^Such type of questions occur, for instance, in the theory of quantum computation, [11,12]. 
Other applications are also possible and will be considered elsewhere. 
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In this case for r = 1, let = while for r > 1, let 




{ (I, J, K) G I for ali p < r and for ali (F, G, H) e 



/6F geG fteH 



The above conjccture in principio can be checked dircctiy using the abovc de- 
scribed recurrence. In reality, for net too large k' s the above recurrence becomes 
impractical. Hence, the conjecture remained unproven for about 36 years since 
its forniulation. A complete solution was found independently by Klyachko [14] 
and Knutson, Tao and Woodward [15, 16] (KTW). Moreover, the infinite dimen- 
sionai generalization of Klyachko results was recently obtained also by Priedland 
[17]. In this work wc shall discuss some results of KTW since, in our opinion, 
they have some physical appeal. In particular, they can be used for solution of 
the following problem. 

Problem 1.3. Is it possible to design a diagrammatic method for de- 
scription of fusion algebra for the Veneziano and Veneziano-like amplitudes ? 
Stated differently, suppose we have the Veneziano (or Veneziano-like) amplitides 
F(A), and F{i'), can we represent the product F{X)F{fi) as 



with coefScients C^^ whose calculation can be complctcly described ? Using 
results by KTW we shall provide an affirmative answer to this question. 

1.3 Organization of the rest of the paper 

In addition to the topics already mentioned we would like to provide a summary 
of the content of this paper. In Section 2 we compare KTW results with those 
obtained much earlier by Heisenberg [18]. Although not present in textbooks on 
quantum mechanics, Heisenberg's originai formulation of quantum mechanics is 
very much in accord with that developed by KTW (who had different purposes 
in mind). For this reason we decided to name the KTW honeycombs ( to be in- 
troduced in Section 2) as Heisenberg honeycombs. The whole discussion of this 
section is motivated by the observation that the honeycomb condition(s) and 
Veneziano condition(s) are mathematically the samc. In Section 3 we review 
some results from Parts II and IH in order to formulate additional problems 
to be discussed in this (expository) paper and later, in its more technical com- 
panion, Ref.[19]. In particular, in this section we formulate problem about the 
fusion rules for the Veneziano amplitudes. Our discussion is not limited to these 
amplitudes however since combinatorially ali scattering processes of high energy 
physics have many things in common. In the samc section we argue (in accord 
with Parts II and III) that combinatoria! considerations alone are sufficient for 




(1.9) 
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recovering the microscopie model leading to Veneziano and Veneziano-like am- 
plitudes. We provide a qualitative evidence that such a model should coincide 
with the standard QCD leaving ali details to [19]. In Section 4 we discuss actual 
calculations of the Littlewood-Richardson (fusion) cocfhcicnts using KTW for- 
malism of honeycombs and puzzles. Although highly nontrivial in their design 
and mathematical justification, these computational tools make calculation of 
fusion coefficients as simplc as possiblc and can bc comparcd in their simplic- 
ity with Feynman's diagrams. The advantage of utilization of such methods 
of calculation becomes especially apparent after we introduce and discuss the 
saturation conjccture and its solution in Section 5. In short, the solution of 
the saturation conjecture allows us to significantly enlarge number of fusion 
coefficients using as an input just few. Moreover, the mathematical proof of 
this conjccture imposcs some uncxpccted extra constraints on fusion cocfhcicnts 
which can be checked experimentally. Section 6 provides necessary background 
for the multiplicative Horn problem to be discussed in Section 7. Section 6 is 
also of independent interest since it provides the most economical and physically 
convincing way to arrive at the classical and quantum Yang-Baxter equations, 
Knizhnik-Zamolodchikov equations, Dunkl operators and Hecke algebra. The 
multiplicative Horn problem discussed in Section 7 cmcrgcs in various physi- 
cal contexts. For instance, it emerges in connection with study of solutions of 
Knizhnik-Zamolodchikov equations, study of spectra of quantum spin chains, 
study of band structure of solids, etc. Although these problems are of no im- 
mediate use for develpment of our string-theoretic formalism, they are logically 
compatible with this formalism and are of independent interest as we explain 
in Section 7. The treatment of these physically interesting problems depcnds 
upon computation of the 3-point Gromov-Witten (G-W) invariants which are 
structure constants of small " quantum" cohomology ring. We enclose the word 
"quantum" in quotation marks since, actually, it should bc callcd "doublé quan- 
tum" or, better, the "deformed quantum" as we shall explain in the text. Com- 
putation of these invariants in Section 7 is non traditional in the sence that it 
is made for people with standard physical education (that is for pcoplc who are 
experts in fields other than string theory). It is hoped, that the provided back- 
ground should be sufficient for proper understanding of current physical and 
mathematical literature on these subjects thus leading to their further uscs in 
condensed matter physics. To keep the main text focused, some computations 
are presented in Appendices A through C. 



2 Heisenberg Honeycomb and Veneziano condi- 
tion 

Following Knutson and Tao [15,16,20] (KT), we would like to describe the 

construction of a honeycomb graph-a prccursor of the Veneziano puzzle-which 
is used for calculation of C^^. In addition, the main mathematical purpose of 
such a honeycomb is to provide constructive solution of the Horn conjecture. 
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In this paper we mainly keep focus on the physical task of calculating C^^. By 
doing so we provide a sketch of how the Horn conjecture was solved with help 
of honeycombs. Details can be found in the Hterature already cited. 

To begin, we need to work out the one dimensionai example (a 1-honeycomb) 
in some detail in order to proceed inductively. It is helpful to replace the 
equation H\ + 'Hfj, = Hi, by the analogous equation H\+H^i+H,y = 0, provided 
that \ + ^ + v = 0. Under such circumstances the inequality z^i < Ai + fj-i 
is replaccd by < Ai + /Xj^ + i^i. Although these results are sufHcient for 
construction of 1-honeycomb, before doing so we would like to put them in 
some historical perspective. 

In his Nobel Prize winning paper [18] Werner Heisenberg made the foUowing 
observation. Influenced by Bohr's ideas, he looked at the famous equation for 
the energy levels difference: 

Lu(n,n- a) = ^(E(n) - E(n~ a)), (2.1) 
h 

where both n and n — a are some integers. He noticed that this definition allows 
him to write the foUowing fundamental composition law 

uj{n — P,n — a — l3) + uj{n, n — 13) — uj{n, n — a — (3), (2-2) 

or, since uj{k,n) — —07(71, fc), the above equation can be rewritten in a more 
symmetric form 

a;(n, m) + w(to, k) + uj{k, n) = (2-3) 

which explains its relevance to the KT results. Most likely, being aware of 
Heisenberg's unpublishcd letter to Kròni^, Dirac in his paper, Ref.[21], of Oc- 
tober 7th 1925 noticed that using the above combinatorial law for the frcquencies 
in the Fourier expansions of observables leads to the multiplication rulc for the 
Fourier amplitudes: 

a{nm)b(rnk) = ah(nk). (2-4) 
By noticing that, in general, 

ah{nk) ^ ba{nk), (2.5) 

he concluded (in accord with Heisenberg's earlier cited letter) that the above 
multiplication rule is characteristic for matrices. Hence, the Fourier amplitudes 
are actually matrices! After this observation, to arrive formally at the famous 
quantization condition 

[x,p]:^ih (2.6) 

using Eq.s (2.4) and (2.5) is quite straightforward. For instance, by using the 
famous Dirac prescription, e.g. see Eq.(ll) of his paper or his book, Ref.[22, 
page 86], 

ih{x,p}p.B. = [x,p], (2.7) 

^Dated by June 5th 1925 (Ref.[18], page 331). 
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where { , }p.b. is the classical Poisson brackct. Such quantization prcscription 
is very formai as Dirac freely admits in his paper, Ref.[21], Section 4. Its cor- 
rectness is based on the historical curiosity which we discuss in detail in Section 
5 in connection with physical (Hciscnbcrg-stylc) solution of the saturation con- 
jecture. In the meantime, we would like to return to Eq.(2.3) in order to analyze 
it from the point of view of modem mathematics. 

Having in mind the KT rcsults [15,16,20], wc would likc to discuss bricfly 
some earlier efforts by Lidskii, Ref. [23], aimed at solution of the Horn conjec- 
ture. Lidskii considered 3n — 1 dimensionai real vector space whose coordinates 
{xi, ...,Xn;yi, ...,yn;zi, ...,Zn) are constrained by the equation 

n 

j:ixj+yj-zj) = 0. (2.8) 

Let Wn be a subspace of R"^" determined by the foUowing inequalities: Xi 
> X2- ■■ > Xn, yi > 2/2 • •• > yn] z\> Z2 - ■■ > Zn- Furthermore, let 5„ be a subset 
of Wn determined by the additional conditions. If a; = (ai, a„; òi, ò„ 
; Ci,...,c„) and x & Sn, then this is possible if and only if there are linear 
Hermitian operators A, B and C acting in n— dimensionai complex space whose 
eigcnvalucs are (ai, a„), (6i, ò„) and (ci, c„) rcspcctivcly, providcd that, 
in addition, A + B = C. The problem lies in describing Sn- Lidskii proves 
that Sn is described by the inequalities given by Eq.(1.7) supplemented by 
the ordering rcquircmcnts on Xi,yi and Zi and the Horn conditions, Eq.(1.8). 
According to Zelevinsky [24], Lidskii actually have not supplied a complete 
proof. Nevertheless, his results apparently have a stimulating effect on the KT 
work. 

Knutson and Tao solved this problem differently in ingenious geometrie way 
by replacing 3n — 1 space with specially chosen 2 dimensionai piane in defined 
by 

R|:=o := {{^^ y, z} e R^ : X + ?y + z = 0}. (2.9) 
The rationalc for this piane can be understood using niethod of induction. In 
one dimensionai case the condition A + /i + = corresponds to some particu- 
lar point in R|-,^q (which wc shall temporarily denote as (A,/i,i^) ). In thcir 
originai work, Ref. [15], KT introduce both the honeycomb tinkertoys and hon- 
eycombs. The tinkertoy is some abstract directed graph placed in R|2=o pl^^^^ 
whose configuration is determined by (cncodcd by) the honeycomb. In ordcr to 
illustrate these concepts we need to establish several additional rules defining 
honeycombs. For instance, with each tinkertoy in R|-^o piane we associate 
another 2 dimensionai (honeycomb) piane (2-plane for short) where lines can 
be drawn only in 3 possible directions: northeast-southwest (Ne-Sw), north- 
south (N-S) and northwest-southeast (Nw-Se) directions. On such 2-plane we 
can place a Y-shaped tripod made of nonoriented labeled semiinfinite edges as 
depicted in Fig.la).The topological Information encoded in such 1-honeycomb 
is used for construction of 1-honeycomb tikertoy, Fig.lb). The location of the 
vertex on such tikertoy R|-;^o pls-nc is determined by the condition X+Ji+P = 
imposed on the vectors A, fi and D in accord with Eq.(2.9). 
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X + /i + = 



Figure 1: The l-honeycomb a) and its 1-honeycomb tinkertoy b) 

Looking at this figure more general rules for n— honeycombs can be devel- 
oped. These are: 

a) there is one-to-one correspondence between the points in ^ piane 

and the vertices in 2-plane; 

b) there should be only finitely many vertices inside the honeycomb diagram; 

c) the semiinfinite lines at the boundary of the honeycomb's boundary are 
allowed to go only in the Ne, Nw and S directions in the 2-plane. 

Remark 2.1. To make sense out of the defining rules for honeycombs, our 
readers are encouraged to look at the following web site: 

http://www.math.ucla.edu/~tao/java/Honeycomb.htmr 
from which they can get an idea of how rules just described are implemented. 

Remark 2.2 By comparing the Veneziano condition, Eq.(1.3b), with the 
definition of R|-, ^ piane, Eq.(2.9), it is clear that the condition, Eq.(1.3b), 

can be brought to the form given in Eq.(2.9). Since the Heisenberg frequency 
condition, Eq.(2.3), is exactly the same as condition in Eq.(2.9) and since it his- 
torically was formulated much earlier, we shall cali just described honeycombs 
as Heisenberg honeycombs. 

Assuming that our readers looked at the suggested web site, we stili would 
like to describe some details about how these honeycombs are constructed in 
order to provide some sketch of KT way of solving the Horn conjecture. 

Although the case of 1-honeycomb is seemingly trivial, already description 
of 2-honeycombs provides much more Information useful for inductive analysis. 
The suggested web link allows our readers to recreate a 2-honeycomb and the 
associated with it 2-honeycomb tinkertoy. Fig.2 provides additional details. 

By discussing these details the rationale for keeping both the honeycomb 
tinkertoys and honeycombs should become apparent. In particular, both the 
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Figure 2: The 2-honeycomb a) and its 2-honeycomb tinkertoy b). The euchdean 
coordinates (x,y,z) of the vertices of such tinkertoy are shown exphcitly 

tinkertoy and honeycomb are uniquely determined by their boundary values, 
i.e. by the prescribed set (A, /x, ly) = (Ai, A2; /i^, /Ì2; ^^i, ^^2) for which the equaUty 
hke that given by Eq.(1.6) is satisfied.Tliat is we must bave 

Al + A2 + +/i2 + i^i + = 0. (2.10) 

In addition, however, for the tinkertoy, the length of the inner edges (determined 
as the Euclidean distance in R^^^q piane between the adjacent vertices) matters 

too. For instance, the Euchdean length of the vector b (may be, up to irrelevant 

factor of \/2) is given by Ai + /Ì2 + i^i > thus providing us with the Ist new 
nontrivial inequality. Clearly, the Euclidean lengths of vectors a and e are 
obtained in the same way and are given respectively by A2 + /Ui + z^i > and 
Ai+/ii+i^2 > 0. In view of Eq.(2.10) these inequalities can be restated as triangle 
inequality for the triangle made of sides whose lengths are A2 — A2, /^i — /U2 and 
L'i — 1/2 respectively. Thus, the Horn conjecture in this case is solved completely 
by simple geometrica! means. For n > 2 the n— honeycomb is determined by 
its 3n boundary values (A, /x, u) = (Ai, A„; /Xi, ui, j/„) subject to the 
constraint analogous to Eq.(2.10). In this case, however, the boundary values do 
not determine the n— honeycomb uniquely. Evcn though they do not determine 
the n— honeycomb uniquely, there is only finite number of different honeycombs 
stili. Evidently, each of these more complex honeycombs will look like that 
depicted in Fig.3 

Already in Part I we noticed that the multiparticle Veneziano amplitude 
can be factorized info product of 4-particle amplitudes, e.g. see Eq. (1.3.28). 
Since with each such factorized amplitude we can associate a 1-honcycomb, it is 
only naturai to expect that the higher order honeycombs can be built from an 
assembly of 1-honeycombs. This indeed happens to be the case and provides a 
compelling reason for our use of honeycombs for description of the combinatorics 
of Veneziano (and Veneziano-like) amplitudes. 

To explain how this happens we need to pay attention to the labeling pattern 
, e.g. that for the 5-honeycomb depicted in Fig.3. One notices that numeration 
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Figure 3: A typical 5-honeycomb 



of external indiccs in the Ne-Sw direction goes from the bottom-up while in 
the Nw-Se direction-from the top-down. In the S direction the numeration goes 
from the right to left. This observation is essential for designing honeycombs. 
We would like to illustrate it by designing, say, a 4-honeycomb. For this purpose 
we have to assemble 4 1-honeyconibs in the way dcpicted in Fig.4. 

What is depicted in this figure is stili a precursor of the honeycomb. To 
convert this precursor info the honeycomb we need to resolve each 4-vertex into 
two 3-vertices as depicted in Fig.5. 

Evcrybody familiar with string theory at this point will be able to recognize 
the famous duality property depicted diagrammatically. c.g. read Rcf.[25], page 
265. Unlikc known particle physics case, such a resolution must be perfornied in 
accord with the rules for honeycombs defined earlicr. It is also ovident that the 
existence of such a resolution is primary cause for the higher order honeycombs 
not to be determined only by the boundary values. 

To illustrate these ideas we apply the rules just discussed to Fig.4 (but 
adopted to 2-honeycombs) . In this case to make the duality of Fig.5 compati- 
ble with the rules defining honeycombs we have to give preference to only one 
type of resolution of 4-vertex. This fact explains why 2-honecomb is determined 
uniquely by its boundary data. This also explains why the higher order hon- 
eycombs are not determined by their boundary data uniquely. The honeycomb 
rules, although quite naturai, stili are not too physically illuminating. To bring 
some physics into this discussion requires several steps. These are described 
below. 



11 



Figure 4: A precursor of a typical 4-honeycomb is made of an assembly of 4 
1-honeycombs 




3 Complex Grassmannians and their relation 
to the Veneziano-like amplitudes 

3.1 Designing partition function reproducing Veneziano 
(and Veneziano-like) amplitudes 

In Parts II and III we constructed the partition function reproducing Veneziano 
amplitudes. In this paper and its companion we would likc to provide additional 
details needed for establishing firm links betwcen this partition function and the 
microscopie model, e.g. QCD. In Part II we used the one-to-one correspondence 
between partitions and the Young tableaux in order to arrive at correct result for 
the partition function. These simple argumcnts were reinforced by much deeper 
results taken from the theory of group invariants for pseudo-refìcction groups. 
In this paper we foUow the same philosophy by providing simple qualitative 
arguments prior to more sophisticated ones. 

We would like to recali some facts from Part II since they will be used 
later, in Section 7, when we shall discuss details of computation of the Gromov- 
Witten invariants. For this purpose we choose a square lattice and place on it 
the Young diagram (tableaux) related to some partition A such that X\- N. To 
do so, we choose some fi x rh rectanglcQ so that the Young diagram occupies 
the left part of this rectangle. We choose the upper left vertex of the rectangle 
as the origin of the xy coordinate system whose y axis (south direction) is 
directed downwards and x axis is directed eastwards. Then, the south-east 
boundary of the Young diagram can be interpreted as directed (that is without 
self-intersections) random walk which begins at (0, — m) and ends at (n,0). 
Clearly, such a walk completely determines the diagram. The walk can be 
described by a sequence of O's and l's, say, for the x— step move and 1 for the 
y—step move. We shall use a notation w(A) for such a walk so that llj{X) := {cui 
,...,ùj]\[) where the random "occupation" numbers luì = 0,1 are analogous to 
those used in the Fermi statistics. Evidently wi +u!2 + • • • + ojn ~ The 
totality M of Young diagrams which can be placed into rectangle is in one-to-one 
correspondence with the number of arrangements of O's and l's whose number 
N is ifi + n. The logarithm of the number J\f of possible combinations of O's and 
l's is just the entropy associated with the Fermi statistic (or, equivalently, the 
entropy of mixing for the binary mixture) used in physics literature. The number 
Af is given hy Af = {rh + njl/rhlnl. It can be represented in two equivalent ways 

(n + l)(n + 2) ■ ■ ■ (h + m) _ f n + m\ 

(m + n)\/m\n\ — ; — \ 

m\ \ m J 

{ih + l)(m + 2) • • • (fi + ?7i) _ / m + n 
ni ~ \ n 

In Part I, Eq-s (I.1.21)-(I. 1.23) explain how the factor Af is entering the 
Veneziano amplitude. Additional significance of this number in connection with 
Veneziano amplitudes is discussed at length in both in Parts II and III. 

^The parameters n and rh will be specificd below. 
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Let now p{N; k, m) be the number of partitions of N into < k non negative 
parts, each not larger than m. Consider the generating function of the foUowing 
type: 



^(fc,m|g)= ^p(iV;fc,m)<?^, 



(3.2) 



N=0 



where the upper hmit S will be determined shortly below. It is shown in Refs.[3- 
5,9] that J^{k, m\ q) k + rh 

k 



pression 



k 



where, for instance. 



k + rh 
rh 



k + m 
fh 



E The 



m 



is a g— analog of the binomial cocfficicnt 



fc + m 
rh 



ex- 



In 



the hterature [3-5,9] this q— analog is known as the Gaussian coefficient. Ex- 
phcitly. 



k 
m 



{q" - - 1) • • • (g* 



-1) 



(3.3) 



From this definition it should be intuitively clear that the sum defining gener- 
ating function jF(fc,TO | q) in Eq.(3.2) should havc only finite number of terms. 
Eq.(3.3) allows easy determination of the upper limit S in the sum, Eq.(3.2). It 
is given by km. This is just the area of fc x m rectangle. Evidently, in view of 

N -k. Using this fact, Eq.(3.2) 



the definition of p{N; k, rn), the number rh 
can be rewritten as: T{N, N — k \ q) ^ 



N 
k 



This expression happens to 



be the Poincaré' polynomial for the complex Grassmannian Gr{rh, k). This can 
be found on page 292 of the famous book by Bott and Tu, Ref . [26j^. From this 
point of view the numerical coefficients, i.e. p{N; k,rh), in the q expansion of 
Eq.(3.2) should be interpreted as Betti numbers of this Grassmannian. They 
can be determined recursively using the foUowing property of the Gaussian co- 
efficients [9, page 26] 



71+1 

k + 1 



n 
k + 1 



k 
m 



(3.4) 



and taking into account that 



n 




= 1. To connect this result with partition 



function reproducing Veneziano (and Veneziano- like) amplitudes we notice that, 



°0n page 15 of the book by Stanley, Ref. [9], 
N{n,k) in positive integers to yi + ... + = i 

number of solutions in nonnegative integers to xi 



one can find that the number of sohitions 
+ fc - 1 
k-1 
n + k 
k 



k is given by 



n is 



whilc the 



.Careful reading 



of Page 15 indicates however that the last number refers to solution in nonnegative integers 
of the equation xq + ... + X/. = n. We have used this fact in Part I, e.g. see Eq. (1.1.21). 

^To make a comparison it is sufficient to replace parameters t^ and n in Bott and Tu book 
by q and A'^. 
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in view of relation m = TV — fc, it is more advantageous for iis to use parameters 
fh and k than N and k. With this in mind we obtairF°l. 

k + m 
k 

L Jq 

= n^^^-^(^'"^N)- (3-5) 

i=l 

This result is of centrai importance since it represents the partition function 
capable of reproducing the Veneziano and Veneziano-hke amphtudes as we have 
explained at length in Parts II and III of our work. In the limit : g — > 1, Eq.(3.2) 
reduces to the number Af as required. To use this information in the context of 
honeycombs we need to remind our readers some basic facts about the Schubert 
calculus 



3.2 Representation theory, Grassmannians, Schubert cal- 
culus and physics of orthogonal polynomials 

The irreducible polynomial representations of general linear group GL{k, C) 
are parametrized by the integer partitions A with at most k parts. Given any 
two such polynomial representations and one can construct the tensor 
product (g) which is expected to be decomposable according to the rule 

V^(E)V'' Gx^V" (3.6) 

into irreducible representations V of GL{k,C). Evidently, since in Eq.(1.9) 
the combinatorics is the same as in Eq.(3.6), the coefhcients C^^ (known in 
literature as the Littlewood-Richardson (L-R) coefficients) should also be the 
same. 

Lct |A| = Al + ... + Afe thcn, it is known [27] that the sum in the r.h.s. of 
Eq.(3.6) is over partitions v for which |A| + |/^| = . Consider the standard 
flag T of complex subspaces C C • • • C C^, where N is related to k 

via k — N — m as before. The Schubert celi ^\ is made out of fc— dimensionai 
subspaces V C with prescribed dimensions of intersections with elenients of 
T . More accurately, they can be described as subspaces for which dim(FnC') = 
wjv +^N-i + • • • +i^N-i+i for i — 1,2, m. Since the notion of the closure Ù\ 
for the Schubert celi is a bit technical (e.g. read page 122 of Ref. [28]) we shall 
skip it without much damage to physics Such closures are called Schubert 
varieties. It happens, that the fundamental cohomology classes a\ = [0,x] 

^"From now on we shall drop the tildas for k and m. 

-'^-'^Effectively, the closure means that the space V defined by dim(V n C*) = uiff + i.^.'jv— i + 
• ■ ■ + aijv— i+i should be replaced by an assembly of spaces for which dim(y n C) > ojjv + 
"^AT-l + ■ • ■ + iUN-i+1 ■ 
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of such varictics forni a Z-basis of the cohomology ring H* {Gr{m, k)) of the 
complcx Grassmannian Gr{m.k) [28-30]. The dimension of such a ring was 

f N \ 

dctcrniincd in the prcvious subscction as , and the multiphcation rule 



k 

for the product of two cohomology classcs givcn by 

o-A-cr^ = EC'a/ìO'>^> (3-7) 

provided that |A| + |/Lt| = \v\ . 

The obtained correspondence can be explained based on "physically intu- 
itive" arguments. Indeed, if the Veneziano (and Veneziano-like) amplitudes are 
perioda of the Fermat hypersurfaces (varieties), as explained in Part I, they 
should be naturally associated with the differential forms generating the coho- 
mology ring for these varieties. When these varieties are embedded into complex 
projective space, where the complex Grassmannian is also embedded (via the 
Pliicker embedding as explained in Part II), the cohomology ring for both of 
these varieties become interrelated. Mathematical details supporting such non- 
rigorous "physical" arguments can be found in the paper by Tamvakis, Ref.[31]. 

The results prescnted above are pretty standard. We would like now to inject 
some physics into them.To do so, we notice that combinatorially the fusion rule, 
Eq.(3.7), is described with help of the Schur polynomials sa(x) [32], that is 
(omitting the a;— dependence) by 

sx-s^ = Y.Cl^s,. (3.8) 

These functions are orthogonal polynomials. That is for partitions and A and 
for the appropriately defined scalar product <, > we obtain: 

(sa, s^) = à\,^J^. (3.9) 

By combining Eq.s(3.8),(3.9) we obtain as well 

<^Am = (sa-s^,s.), (3.10a) 

which, in view of Eq.(3.7), is equivalent to 

CAM = (^A-CT^,a.). (3.10b) 

Although this result is very important from the point of view of algebraic ge- 
ometry (since it describes intersection of Schubert cycles), to use it physically 
requires more work as we would like to explain now. 

In Sections 7 and 8 of Part II we discussed why the Schur polynomials 
associated with KdV hierarchy (and, hence, with the Virasoro algebra (through 
method of coadjoint orbits) are not relevant for Veneziano amplitudes. At this 
point we are ready to provide additional explanation why this is so. 

First, we would like to recali the defìnition of one of the basic symmetric 
function m\ [32]. For this we have to define the monomials, e.g. x"^ = x\^x'^ ■ ••, 
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associateci with partition A. With thus defincd monomials, m\ is just the sum of 
these monomials made of ali possible permutations oi x' s. When this definition 
is combined with the results of Part I desribing Veneziano amplitudes as period 
of Format varictics, it should bccome elcar that the fuUy symmetrized Veneziano 
amplitude can be obtained by using rnx in the numerator of the period integrai. 

Remark.3.1. It is known that m\ are eigenfunctions of the Calogero- 
Sutherland (C-S) model [33]. It is also known that 2 dimensionai QCD is re- 
ducible to the C-S model [34,35]. Thus, it should be not too surprising that 
Veneziano amplitudes had been successful in describing scattering of mesons. 

Second, since wc are intcrcstcd not only in Veneziano amplitudes but in 
general combinatorial properties of the scattering amplitudes of high energy 
physics, we would like to develop things a bit further. For instance, in the next 
subsection wc shall demonstrate that combinatorial data contained in Veneziano 
amplitudes are quite sufìicient for calculation of Cj^^. From this fact it is easy 
to make a mistake and to use the Schur polynomials sa(x) in the subsequent de- 
velopments. This is known and well dcveloped pathway to the traditional string 
theory formalism. One should keep in mind however that, like in ordinary quan- 
tum mechanics, any symmetric function can be Fourier decomposed into Fourier 
series whose basis is madc of orthogonal polynomials. In view of Eq.(3.9), the 
Schur polynomials provide a basis for such an expansion but this basis is not 
unique. There are other orthogonal polynomials which can be used for such 
a purpose [36,37]. Very much like in quantum mechanics, where ali exactly 
solvable problems possess a complete orthogonal set of eigenfunctions, differ- 
ent for different problems, one can think of the corresponding exactly solvable 
(many-body) problems associated with orthogonal polynomials, also different 
for different problems. More interesting, however, is to be able to solve the 
inverse problem. 

Problem 3.2. For a given set of orthogonal polynomials find the corre- 
sponding many-body operator for which such a set of orthogonal polynomials 
forms a complete set of eigenfunctions. 

In view of the Remark 3.1., this task can be solved completely in the case 
of Veneziano amplitudes. In this paper and its companion we make an attempt 
at providing more general outlook at solution of the Problem 3.2. We hope, 
that by rasing these issues more works will follow enabling to solve this problem 
completely. 

4 Designing and solving Veneziano puzzles 

4.1 General remarks 

With the background just provided we are ready to connect the combinatorics of 
honeycombs with computation of the L-R coefBcients C^^. As a by product we 
shall introduce another honeycomb-related construction for calculating these 
coefficients which KTW cali a "puzzle" [16]. Combinatorics of honeycombs is 
connected with the L-R coefficients in view of the foUowing theorem, Ref.[15], 
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if' t 



Figure 6: A typical example of a graphical calculation of the L-R coefBcients 
with help of 3-honeycombs 

page 1053. 

Theorem 4.1. Lei X, fi and v he three pre assigned (boundary) partitions 
(e.g. like those depicted in Fig.3 for the 5-honeycomb) for the k~honeycomb. 
Then the number of dijferent honeycombs with such pre assigned boundary con- 
ditions is given by the L-R coejjìcient C^^ 

Although the above theorem hints at physical relevance of honeycombs, e.g. 
for practical calculation of the L-R cocfficients, actual use of honeycombs for 
such a purpose based on the Information provided is somewhat problematic. We 
would like to correct this deficiency now. Firstly, foUowing KT [15], we would 
like to illustrate general principles by using a simple example. In particular, in 
the case of 3-honeycomb the decomposition of the tensor product 

^(24,0) ® ^(2,1,0) = ^(4,2,0) © ^(3,2,1) ® ^(4,1-1) ® ^(3,3,0) © ^(2,2,2) 

is graphically depicted in Fig.6. 

Evidently, the L-R coefìicients can be read off from such a decomposition 
straightforwardly. Moreover, as Fig.6 indicates, in actual calculations of these 
coefficients it is sufìicient to use the precursor, e.g. see Fig.4, rather than the full 
blown 71— honeycomb. This observation is especially helpful for the Veneziano 
(and Veneziano-like) amplitudes in view of already noticed factorization prop- 
erty provided by Eq. (1.3.28). It is more questionable if we are interested in the 
most general form of multiparticle scattering amplitude compatiblc with the 
energy-momentum conservation laws. 

For the sake of such generality, we would like to discuss yet another method 
of computation of the L-R coefììcientJ^. This method requires designing and 
solving puzzles associated with honeycombs. The simplest puzzle associated 
with 3-honeycomb is depicted in Fig.7. 

^^There are many methods of calculating these coefficients. In fact, the number of these 
methods is so large that even the most rcputable monographs, e.g. Ref. [37], are unable to 
provide the complete hst. We only mention K-T hives to be discussed in Section 5.3. which 
are just a shght modifications of the Berenstcin-Zclevinskii (BZ) patterns nicely described in 
Ref. [37], page 437. We discuss KT variant of constructing the L-R coefficients mainly because 
of the factorization property, Eq.(I.3.28), of the Veneziano amplitudes. 
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4 



Figure 7: The labeling order of vertices of the 3-honeyconib tinkertoy a) is in 
one-to onc corrcspondcncc with labehng of triangles (precursors of) in a puzzle 
b) associated with such a tinkertoy 



This picture only provides a hint that honeycombs and puzzles are intercon- 
nected but not much insight into rationale for switching from honeycombs to 
puzzles. We would like to discuss this rationale now. For this purpose, we need 
to use the correspondence between the partitions and directed random walks 
already discussed. For each partition triple A, fj, and i/ there is its realization in 
terms of such walks: w(A), and u){i'). Consider a particular directed random 
walk. It is described by the Fermi- type variablc loì{X) such that the constraint 
^^■^ Wi(A) = m = N — k holds. Consider now an cquilatcral triangle whose 
sides are divided into N segments of equal length. Furthcrmore, let us put 
these segments in correspondence with llIì{X),llIì{p,) and llì-ìÌv), respectively for 
each side of the triangle. Finally, consider the set of puzzle pieces depicted in 



They are made of cquilatcral triangles and rhombi whose sides ali have the 
same lengths cqual to that of the segmcnt on the side of the larger triangle. 
Since thcsc puzzle picccs are labeled, the task is to fili in the large cquilatcral 
triangle with the puzzle picccs providcd that thcsc picccs can bc rotatcd but 
not rcflcctcd whcn they are uscd to solve the puzzle. The final rcsiilt looks like 
that depicted in Fig.9. 

By analogy with Theorem 4.1. it is possible to prove the following theorem. 

Theorem 4.2. (Knutson-Tao-Woodward [16]) Let partitions A,/i and v he 
encoded by random walks a;(A).u;(/i) a,nd lo{v) whose particular realization is 
described in terms of 's and l 's on the sides of the cquilatcral triangle encoded 



AVA y' 




Figure 8: Basic building blocks of a puzzle 



Fig.8. 
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Figure 9: A typical assemblee! puzzle 

from left to righi, on each of its sides in a clockwise arder. Then, the n,um,her of 
puzzles constructed with such boundary data equals to the L-R coefficient C^^. 

The easiest way to prove this thcorcm is through graphical bijection between 
the honeycombs and puzzles which we would like now to describe. At this point, 
in view of Fig.7, we know already that such a bijection does exist. Using results 
of KTW [16], it can be made more accurate now. For this purpose, the following 
steps should be made: 

a) one should place a solved puzzle on the R|-, ^ piane in such a way that 

the bottom right corner is at the origin. Next one should rotate this puzzle 
around origin by 30*^ countcrclockwisc; 

b) at each boundary segment, which is labeled by l(l-region), attach a rhom- 
bus (outside the puzzle), then another (parallel to fist) and so on ad infimum. 
Fili in the rest of the piane with 0-trianges (see Fig.lO): 

c) deflate thus obtained extended puzzle, while keeping the right corner at 
the origin. The result will be the honeycomb whose vertices originate from such 
deflated 1-regions and whose edges are labeled by the thickness of the originai 
rhombus region. 

Fig.lO illustrates such described reduction procedure. 

This complctes our description of KTW puzzles and the associated with 
them Heisenberg honeycombs. 

4.2 Some comments about solution of the Horn conjecture 



The results described previously are incomplete without further discussion 
of the Horn and saturation conjectures. In this subsection we would like to 
discuss some additional details related to the Horn conjecture. These may be 



20 



Figure 10: Successive stages of conversion of a puzzle into honeycomb 



Figure 11: Designing largar honeycombs (honeycomb tinkertoys) from the 
smaller ones. The origin (0,0,0) is marked for each case by the black dot 

of some importance in organizing experimental data for the mass spectrum of 
hadrons. 

Since the case of 2-honeycombs can be solved completely, the foUowing prob- 
lem emerges. 

Problem 4.3. To what extent can one use methods developed for 2- 

honeycombs to obtain similar typc incqualities for more complex honeycombs? 

We can develop our intuition by fooking at Fig.4 and asking a question: is it 
permissible to make, say, 3-honeycomb out of 2 and 1-honeycombs as depicted 
in Fig.ll. 

Stated more formally, suppose we have h and h' honeycombs what is the 
meaning of their direct sum h ® h'I li we are talking about n— honeycomb h 
and 777— honeycomb h' , thcn the direct sum h ® h' must evidently correspond 
to the direct sum of n x n and m x m matrices combined together to form an 
(m + n) X (m + n) block-diagonal matrix. With such clarifications, it is possible, 
foUowing KT, to give a purely geometrie proof (involving n-honeycomb) of the 
inequality Aj + jij -\- > (known already to Weyl, Ref.[10]) provided that 
i + j + k = n + 2. For n = 1 this inequahty is obviously correct. For n > 1, 
the requirement i + j + k = n + 2 should hold while in order to prove that 
+ A*j + i^fe > 0, the "physical" arguments can be used as foUows. Consider 
some vertex P inside the n— honeycomb as depicted in Fig.l2. 
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Figure 12: Using honeycomb for proving Weyl's inequality 



Next, we connect this vertex with the boundary lines marked respectively 
as Ai , iij and Vk ■ Assume that our honeycomb is made of wire and assume that 
Constant currents with intcnsity A^, /i^- and Uk are apphcd at the boundary hnes 
whilo the rcst of the boundary hnes serve as sinks (i.e. they are grounded). 
Remembering the Kirkhoff's rule for each vertex (that the total algebraic sum 
of ali currents cntering a givcn vertex must be zero) we apply this rule to the 
selected vertex P. If now the edges emanating from P are labeled respectively 
by A-, and then, the Kirkhoff rule requires: A - + i^j + =0. At the same 
time, since the current goes to other vertices too, it should be clear that: Aj > A^, 
A*j ^ Mj- ^-nd i^fe > i^Jj. From here we obtain: Aj + fij + > 0, in accord with 
Weyl. Since we would like to stay focused on our immediate (physical) tasks, we 
refer our readers to the already cited literature for additional details on solution 
of the Horn problem. 

5 Solution of the saturation conjecture: from 
Heisenberg to Knutson and Tao and beyond 

5.1 Statement of the problem 

Before providing exact definitions, following KT we would like to discuss some 
general issues for the sake of presenting them subsequently in a proper physi- 
cal context. In particular, if there exist Hermitian matrices A,B and C with 

eigenvalue sets A, fi, and u respectively, one can associate to the matrix equation 
A + B + C = the symbolic relation of the type 

Affl/ifflz/~cO<s=^ A + /i + z/ = 0. (5.1a) 
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Here the subscript c mcans " classical" . Although the hermiticity of matrices 
A,B and C makes them suitable for quantum mechanical interpretation, KT 
introduce another relation 

Affl/xffli^~, A + A* + J/ = (5.1b) 

where the subscript q means "quantum". KT consider "classical" problema 
as those involving Hermitian matrices with real spectrum while they consider 
problems as "quantum" if the spectrum involves only integrai A,/i, and u's. 
Based on these definitions, KT formulate and solve the following theorem. 

Theorem 5.1. (Knutson-Tao [20]) Let A,/i, and v he weakly decreasing 
sequences of n integers. Then, if Eq.(5.1b) holds, Eq.(5.1.a) holds as well. On 
another hand, if Eq.(5.1a) holds, then, there exists an integer N > such that 
N\ ffl Nfi ffl Ni' 4=^> A ffl ^ ffl 0. (That is asymptotically, Le. for some 
large N's, quantum and classical results coincide) 

Based on this result KT formulate and prove the following 

Conjecture 5.2. (Saturation conjecture) The above classical-quantum 
equivalence persist even for A'' = 1. 

Remark 5.3. The above conjecture can be restated in terms of the L-R 
coefRcients C^^. Following Fulton, Ref.[38], page 238, we expect that if A, fi, and 
v are a triple of partitions and C^^^^ ^ 0, then C^^ ^ as well. In particular, 
if C^^ = 1 we should expect C^^j^^ = 1. 

It should be noted, however, that in the case if C^^ ^ 1 there is no rea- 
son to expect that C^'^p^^ = C^^. Since such an observation can be checked 
experimentally, it makes sense to discuss it in some detail in this section. Math- 
ematically, the proof of saturation conjecture was made not only by KT but by 
several other authors as well . They are listed on page 238 of Ref . [38] . In this 
section we would like to discuss the physics of this proven conjecture based on 
arguments used by Heisenberg in his key paper on quantum mechanics, Ref. [18]. 
Some auxiliary results are presented in Appendix A. 

5.2 Heisenberg's proof of the saturation conjecture 

We begin by noticing that use of Hermitian operators in quantum mechanics 
is motivated by the requirements that the observables which these operators 
represent are real mmibers. In particular, for isolated stable physical system 
the spectrum of its eigcnvalues should be real. This fact is in apparent contra- 
diction with the KT definition of what is "classical" and what is "quantum". 
For instance, the famous Hydrogen atom energy spectrum is known to behave 
as En ~ l/n'^.n = 1,2,... The contradiction, nevertheless is only apparent. 
Before going into detailed explanations, we recali that famous semiclassical ap- 
proximation of quantum mechanics relates quantum results to classical in the 
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limit of large quantum numbers. This fact can be taken as physical proof of 
Theorem 5.1. This makes sense only if we accept that the Hermitian operators 
producing real spectra bave something in common with the "classical" world. 
Alternatively, we can try to prove that such a spectra cannot belong to any 
quantum mechanical system. Clearly, such a proof will stili be insufhcient to 
place such type of spectra into "classical" world since in classical world there 
are no operators and everything commutes. Hence, we would like to approach 
the saturation conjecture somewhat pragmatically using physical arguments. 

For this purpose, we would like to bring to attcntion of our readers some 
important quotations from the classical hook by Dirac, Ref.[22]. On page 177 of 
this book one reads: " In fact it was the idea of replacing classical Fourier com- 
ponents by matrix elementJ^^lwhich lead Heisenberg to the discovery of quantum 
mechanics in 1925. Heisenberg assumed that the formulas describing the inter- 
action with radiation of a system in the quantum theory can be obtained from 
the classical formulas by substituting for the Fourier cornponents of the total 
electric displacemenlF^ of the system the correspondinJn matrix elements" . 
Further in the text Dirac elaborates on this quantum-classical correspondence. 
Specifically, on pages 245-246 we read: " Thus, the elementary theorj0...in 
which the radiation is treated as an external perturbation[3 gives the correct 
vaine for the absorption coeffìcienlF^. This agreement between the elementary 
theory and the present theory could be inferred from general arguments. The 
two theories differ only in that the field quantities ali commute with one an- 
other in the elementary theory and satisfy definite commutation relations in 
the present theor^P^. and this difference becomes unimportant for strong fields. 
Thus the two theories must give the same absorption and emission when strong 
fields are concerned. Since both theories give the rate of absorption proportional 
to the intensity of the incident beam, the agreement must hold also for the weak 
fields in the case of adsorption. In the same way the stimulated part of emission 
in the present theory must agree with the emission in the elementary theory" . 

We brought such extensive quotations from Dirac only to emphasize that, 
actually, at least in some cases, the quantum-classical correspondence can be 
pushed way down into the seemingly quantum domain thus providing a "proof 
of the saturation conjecture. Since the existing literature on quantum mechan- 
ics (including Dirac's book) for some reason does not discuss these issues in 
sufficient detail, we would like to provide such details in this section. 

Following the logie of Heisenberg's originai paper, Ref.[18], Eq.(2.1), when 
it is treated classically, can be rewritten as 

Lo[n^a) — auin) ~ a——— . (2. la) 

h on 



"E.g. see Eq.s (2.4), (2.5). 
^^The dipole moment. 
^^I. e. quantum. 

^^That is completely classical (e.g. see Appendix A) 
^^That is classically! 

^^Which is calculated quantum mochanically 
-^^That is quantum. 
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Next, Heisenberg notices that, actually, the famous Bohr-Sommerfeld (B-S) 
quantization rule 

pdq = nh, ri = 0, 1, 2, (5-2) 

is not exact ! It is determined only with accuracy up to some Constant (unknown 
at the time of his writing). He argues, that if such a Constant would be known, 
this rule would become exact, that is valid for any n's. From the point of view 
of our present understanding of quantum mechanics his intuition was correct: 
the old fashioned Bohr-Sommerfeld rule is valid rigorously in the limit of large 
n's while the calculation of the Constant can be done, for instance, with help 
of either the WKB or of much more sophisticated theory of Maslov index 
[39]. These arguments although plausiblc are superficial ncvertheless as can 
be found from the hook by Arnold, Ref.[40], page 246. From it wc find that 
already at the classical level the adiabatic invariant § pdq is determined only with 
accuracy up to some Constant. This observation makes Heisenbeg's arguments 
less convincing. 

In particular, foUowing Heisenberg we claim that if the B-S quanization 
rule, when corrected, makes sense fully quantum mechanically, one can get, in 
principle, the additional Information out of it. For this purpose Heisenberg 
introduces the Fourier decomposition of the generalized coordinate i.e. 



oo 



q{n,t) — aoi{n) eyi'p{iu){n,a)t) (5-3) 

a— — oo 

where we used Eq.(2.1aj^ thus causing us to keep our calculations with respect 
to the some pre assigned energy level n (e.g. see Appendix A). The velocity can 
be readily obtained now as 

oo 

q{n,t) = iaa{n)Lu{n,a) exp{iuj{n,a)t) (5.4a) 



a — — oo 



so that calculation of the velocity square averaged over the total period is given 

by 

/oo 
[q{n, t)f dt ^2Tr ^ |a„(n)|^ a)^ (5.4b) 



Q — — oo 



At this point it should be noted that the originai of Heisenberg's paper, Ref. [18], 
contains (perhaps) a typographical error: instead of having uj(n, a)^ Heisenberg 
writes U!{n,a). This fact was noticed by the editors of his coUected papers, 
Ref. [18]. Now we use this result in the B-S quantization rule, i.e. we have 

j) pdq — ^ mqdq = j) mq^ dt = 2T:m \aai'n)\'^ uj{n, a)'^ = nh + const. 

(5.5) 



Q — — oo 



^"In the originai Heisenberg keeps aiL>{n) instead of u!(n,a). Our notations happen to be 
more convenient as we shall demostrate shortly. 
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Ncxt, Heisenberg argues as follows. Since the const is unknown, it is of interest 
to obtain results which are constant-independent. At the same time, since the 
result, Eq.(5.5), is assumed to be exact, we have to use instead of scalars |aa(n)|^ 
the matrices in accord with Eq.(2.4). This would Icad us to matrices of the type 
\a{n,n + a)\^ and \a{n,n — a)\^ depending on the actual sign of a. In addition, 
he silently had assumed that the n— dependence of the amphtudes is much 
wcaker than that for the frequcncies Cij(n, n + a) and Lu{n,n — a) so that it can 
neglected completely. Under such conditions he treats n as continuous variable 
and differentiates both sides of Eq.(5.5) with respect to n thus obtaining (recali 
that u){mn) = —w{nm)) : 

oo 

h = Anni ^^{|a(n, n + a)|^ u>{n, n + a) — \a{n, n — a)\^ oj{n, n — a)}. (5.6) 
«=0 

The validity of this result depends upon the additional assumption about the 
ground state energy. If no represents such a state, then one must require that 

a(no,no-a) = Va > 0. (5.7) 

In Ref.[18] Heisenberg acknowledges that this result was inspired by the 
earlier result of Kramers who calculated the induced dipole moment of electrons 
in atom assuming rules of quantum mechanics (in fact before it was ofììcially 
inaugurated !), e. g. see Appendix A. Results of the Appendix A then lead us 
directly to the famous commutation rule 

[x,p]=ih. (5.8) 

Heisenberg argues that his reasonings are correct since the frequency of the 
incoming (scattered) light is much higher than that for characteristic "rota- 
tional" frequencies in the atom so that the electrons can be treated as "free" 
and independent. 

The discussion we just presented is aimed to underscore the differences be- 

tween physical reality and mathematical corrcctness. It can be considered as 
the Heisenberg-style proof of the saturation conjecture and provides us with 
rationale for discussion of an alternative formulation of quantum mechanics (to 
be presented in the next section). Before doing so we would like to discuss 
the saturation conjecture and its proof in connection with results of our earlier 
published Parts II and III. This will enable us to make some physical sense out 
of recently obtained mathematically interesting results. 

5.3 Combinatorics of L-R coefRcients and the Ehrhart poly- 
nomial 

In this subsection we would like to address the following problem: suppose we 
are given a L-R coefhcient C^^, can this Information bc used for calculation of 
^nxnlt Although calculations of L-R coefHcients have a rather long history [41], 
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the full answcr to tliis qucstion was obtaincd only quitc rcccntly in connection 
with positive solution of the saturation conjecture. It should be noted that 
although the attempts in this direction were made a bit earlier by Berenstein 
and Zclcvinsky [42], the actual numcrical rcsults wcrc obtaincd much latcr by 
King et al [43]. These authors were inspired by the results of KT, Refs. [15,16], 
where, in addition to the honeycomb model, the hive model was introduced 
which we have not discusscd thus far. An n-hive is a triangular array of numbers 
ttij with < i,j < n. Say, for n = 4 a typical arrangement looks as foUows 

o-oo 

aio doi 

«20 «11 «02 

«30 tl21 «12 «03 

O40 «31 «22 0,13 «04 



Such an n— hive is an integer hive if ali of its entries are non-negative integers. 
The numbers a,j in the hive are subject to the hive conditions given symbolically 

ab b d 

by RI: , ; R2: b c and R3: implying ò + c > a + d for 

ed , a c I- ^ c - 

d 

a, h, a and d. bcing the ncighboring cntrics in the n— hive. Wc shall cali such type 
of inequalities the hive condition (HC). Based on these results, the foUowing 
definition can be made 



Definition 5.4. A L-R hive is an integer hive satisfying the HC for ali 
constituent rhombi R1-R3. For such a hive the border entries are dctcrmined 
by partitions A, and v in such a way that aoo = 0, aoj = Xi + X2 + ... + Xj ,j = 
1, 2, n, aio = Ui + U2 + ... + i^i, ? = 1, 2, n, ak,„-k = aon + Mi + /"2 + ••• + Mfe 
, k = 1,2, n, provided that |A| + = \v\ and the number of parts Z(A), 
and /(!^)in partitions A, n and v is bounded by n. 

Based on this definition, and motivated by Theorems 4.1. and 4.2. Fulton 
[44] proved the foUowing theorem 

Theorem 5.5. The L-R coefficient C^^ is the number of LR hives with 
border labels determined by A, jx and v. 

Alternative (simpler) proof of this theorem as well as many other useful 
results can be found in the recent paper by Pak and Vallejo [45]. In both cases 
the proof is based on careful solution of the rhombus constraints (or HC) for the 
integer hives. In the case of n-hive there are m = (n— l)(n — 2)/2 interior vertex 
labels for ajj.The corresponding set of linear constraints with integer coefficients 
defines a convex rational (not integrai! as in Part III) (hive) polytope V living 
in R™. As in Part III we can define the Ehrhart polynomial AT, for such 
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polytope and the associateci with it generating (partition) function T(V^x) via 

oc 

T{V,x) = rn)x^ . (5.9) 

In the present case, <P(Ar, m) =C^li^^ E3. Unlikc the integrai polytope V for 
which the generating function J^{V,x) can be written in closed forni given by 
Eq. (III. 1.14), in the present case, since the polytope is only rational (that is not 
ali of its vertices lie at the nodes of Z™), such closed forni universal result for 
T{V , x) cannot be used. Using niethod of quivers, Derksen and Weyman [46] 
were able to prove that, nevertheless, the universal forni given by Eq.(III.1.14) 
stili holds (provided that the diniensionality m of the lattice Z™ is replaced by 
m and the indeterminate x is replaced by a;" where both fh and a are some 
known (in principle) nonnegative integers. 

In view of such an interpretation of stretched L-R coefficients, the whole 
chain of arguments of Part III can be used practically unchanged. This makes 
Eq. (III. 4. 30) especially relevant demonstrating that combinatorics of the L-R 
coefficients can be obtained field-theoretically in the spirit of earlier treated 
Witten-Kontsevich model [47] . These results will be extended further in Section 
7 in which we discuss (among other things) the Gromov-Witten invariants. 

Remark 5.6. Obtained connections between the L-R coefficients C^^ and 
their inflated counterparts C^'^jq^ is very important from the experimental point 
of view. Indeed, by looking at Eq.(1.4) and recalling its relevance to Veneziano 
amplitudes described in Parts II and III, one would naively expect that for any 
N one should have C^^ = C'^'^j^^ . Mathematics tells us that this may happen 
only if C^^ = 1. It remains to check experimentally if the constraint C^^ = 1 
is a valid physical constraint. These arguments clearly not restricted to the 
Veneziano-type amplitudes and should be taken info account for any amplitude 
of high energy physics. 

Problem 5.7. In the case if the hypothetical scattering processes requires 
use of supersymmetric QCD, how supersymmetry can be detected through study 
of combinatorics of experimental data, say, of the LR fusion coefficients 



6 From Heisenberg back to Maxwell, Tait and 
Kelvin 

6.1 General remarks 

In view of the results of previous section and those of the Appendix A we would 
like to formulate the foUowing problem 

^^A typical example is shown on page 14 of Ref. [43]. 

^^Even though methods developed in Part III include use of supersymmetry, it was demon- 
strated in Parts II and III that its use is not essential. Hence, one should not confuse these 
results with the supersymmetric extension of the underlying microscopie QCD model and with 
calculations of observables for such type of models. 
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Problem 6.1. Given the experimental origin (discusseci in Appendix A) 
and the associated with it inevitable computational approximations leading to 
the basic commutation rule, Eq.(5.8), is it possible to develop quantum me- 
chanics without this rule being put at its foundation? Alternatively said, can 
we recover this basic rule without use of light scattering experiments and the 
B-S quantization rule? 

Remark 6.2. Since ali experimental data for any quantum mechanical 
system (atom, molecule, solid, etc.) are spectroscopic in nature we know about 
the system as much as the combinatorics of experimental data provides. From 
such point of view there is not much difference between, say, biological problems, 
computer science problems, astrophysical problems, etc. and those in the high 
energy physics. 

A superficial answer to just posed problem can be madc like this: should the 
above rule be wrong we would not be able to recover the spectrum of Hydrogen 
atom with such an amazing accuracy using established methods of quantum me- 
chanics. This remarkable agreement between theory and experiment is possible 
only if such commutation rule is correct. Clearly, it is correct! Nevertheless, 
we can stili argue against its up-front use as follows. Heisenberg uses the B-S 
quanization rule (perhaps adjusted) to obtain his results. This rule makes sense 
only if classically there is a complete separation of variables done with help of 
the Hamilton- Jacobi formalism. When this happens, the system is considered 
to be completely integrable. Hence, any completely integrable system is just 
the set of independent harmonic oscillatori. The Hydrogen atom is surely a 
good candidate for such procedures but what about the Hclium ? The B-S rule 
cannot be applied strictly speaking already to the HcliunQ so that Heisenberg' s 
chain of reasoning leading to the commutator, Eq.(5.8), formally breaks dowro. 
Besides, since the B-S quantization cannot be used for spin quantization (since 
formally there is no classical analog of spin, i.e. B-S rule does not account for 
the half integers), the spin has no place in the Schròdinger's formalism. Since 
Schròdinger have demonstrated the equivalence of his formalism to that devel- 
oped by Heisenberg as described in Dirac's book, Ref. [22], apparently, there is 
no room for the spin in the Heisenberg formalism as well. Surely, this happens 
to be only apparently true as we would like to explain now. In doing so we do 
not need to use the relativistic formalism developed by Dirac. 

^''This explains why Heisenberg was able to do his formai differentiation (over n) of Eq.(5.5) 
and arrived at correct result. This also explains why KT cali system "quantum" if it has an 
integrai spectrum according to the B-S quantization prescription. 

^^A very interesting detailed discussion of this fact is given in the monograph by Max Born, 
Ref. [48], pages 286-299, published in 1924, i.e.prior to the officiai birth of modem quantum 
mechanics. 

^^In the paper by Pauli and Born written in 1922 [49] it is noted that Bohr conceded 
that only the Hydrogen atom is quantizable but the rest of atoms are not. Therefore the 
spectral lines of elements other than Hydrogen must be noticeably wider. This expectation 
is in disagreement with what is observed specrosopically. The spectroscopic data for most of 
elements of periodic table were available already in 1905 [50], that is long before the quantum 
mechanics was formulated. 
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6.2 Heisenberg's paper revisited 

We begin with discussion of some consequences of Heisenberg's results foUowing 
thel925 paper by Dirac, Ref. [21]. We selected this paper in view of its remark- 
able completeness: ali quantum mechanical formalism used today can be traced 
back to this papeiF^. Dirac acknowledges, though, that bis paper was written as 
consequence of Heisenberg's results. In particular, the famous Dirac quantiza- 
tion rule, Eq.(2.7), is just restatement of the results by Heisenberg. As good as 
it is, its use is questionable in general. Indeed, in comments to bis Eq.(ll) Dirac 
States that the difference of Heisenberg's products of two quantum observables 
X and y is equal to the (classical!) Poisson bracket of their classical counterparts 
multiplied by ^ or, symbolically, 

xy - yx ^ ^{x,y}p,b, (6.1) 

This expression makes sense for a; ^ i; and y ^ p . But, in general, for 
arbitrary classical observables x and y we are dealing with the Lie algebra which 
requires this Poisson bracket to be expanded into linear combinations of classical 
obseravbles so that the l.h.s and the r.h.s of Eq.(6.1) do not match. Hence, we 
come back to the Heisenberg-Kramers result presented in the Appendix A as the 
only justification. This difRculty was recently noticed and discussed in the book 
by Adler who suggests to treat classical Poisson bracket quantum mcchanically 
in the style of Heisenberg, e.g. see Eq.(1.13b) of Ref.[52]. 

In this paper, we choose another way to by pass the Heisenberg-Dirac quan- 
tization prescription. For this purpose, we would like to make few additional 
comments regarding traditional formulations. FoUowing Dirac we introduce the 
evolution operator U (t) bringing the initial state wave function ipQ to its final 
state ip{t), i.e. ip{t) = U{t)ipQ. For the time-independent Hamiltonian H the 
formai solution of the Schròdinger-type equation 



dt 



ih^U-\t) ^ HU-\t) (6.2) 



is known to be given by Ù^^{t) = cxv(—j-Ht)W\ Now, Heisenberg considered 
the quantum Fourier transform by replacing the usuai Fourier amplitudes by 
matrices, i.e. he used quantities like a(rnn)exp{jrLu(rnn)t) . In the modem lan- 
guage this can be rewritten as foUows. Let O be some quantum mechanical 
operator whose evolution is described by U{t)OÙ^^{t) — 0{t). This operator 
leads to the matrix elements : < m \ O \ n > exp{jruj{mn)t) with uj{mn) de- 
fined by Eq.(2.1) with < m | and | n > being time-independent wave functions 
of the Hamiltonian H. Clearly, if the observable 0{t) is an identity element 



^^It should be noted, nevertheless, that Dirac's paper [21] was received by the Editoria! 
office on 7th of November of 1925 while on November 16th of the same year the paper by 
Borri, Heisenberg and Jordan, Ref. [51], was registered by the Editors. It contained practically 
the same results as Dirac's paper and many other results in addition. 

^^We write instead of U to be in accord with mathematical literature. This will be of 
immediate use shortly below. 
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in the algebra of observables, we obtain : < m{t) \ n{t) >—< m \ n >, where 
I n{t) >= U{t) \ n > . The requirement for the observables to be real leads 
to the Hermitian type operators whose eigenfunctions are mutually orthogonal. 
Under such conditions we may or may not require these mutually orthogonal 
functions to be normalized to 1. By doing so we do not insist on the proba- 
bilistic interpretation of quantum mechanics. Such an interpretation emerges 
anyway within quantum statistical mechanics. 

In view of earlier posed Problem 6.1. it makes sense to replace the Dirac 
quantization rule, Eq.(6.1), by the requirement of orthogonality for the wave 
functions. Under such a rule we need to have a supply of orthogonal functions 
(if the spectrum is countably infinite) or orthogonal basis in some complex fi- 
nite dimensionai vector space. In the case of orthogonal functions, it is known 
that ali one variable orthogonal functions used in quantum mechanical exactly 
solvable problems are obtainable from the one variable Gauss-type hypergeo- 
metric functions [53]. These functions are expressible in the form of period 
integrals. The Veneziano and Veneziano-like scattering amplitudes considered 
in Part I belong to the same family of hypergeometric- type period integrals 
initially considered by Aomoto [54] and subsequently by many others [55] . The 
cohomological meaning of such integrals is explained in detail in Ref. [56] . By 
the principle of complementarity ali many-body exactly solvable quantum me- 
chanical problems should be related to the hypergeometric functions of multiple 
arguments. More importantly for us is that these hypergeometric functions pro- 
duce sets of ali known orthogonal polynomials replacing one- variable orthogonal 
functions of usuai quantum mechanicj^. Hence, they are also period integrals. 
A nice summary of developments in this area can be found in Refs. [57,58] . The 
finite dimensionai cases (including spin) technically present no difhculties under 
such circumstances. 

At this stage we are ready to provide additional arguments in support of 
our point of view on quantum mechanics. These arguments will be also of use 
in the next section. Traditionally, in the Heisenberg interpretation of quantum 
mechanics equations of motion for the operators Oi (t) can be obtained by simple 
differentiation of U{t)OiU~^{t) = Oi{t). This procedure formally leads to the 
Heisenberg's equation of motion 

dò - - 

ih^^[0,H] (6.3) 

for the operator O. The rationale for such writing comes from the analogy of 
this equation with that known in classical Hamiltonian mechanics. This makes 
sense only if the Dirac quantization prescription makes sense. But it does not as 
we just discussed! Instead of repairing this situation using known mathematical 
methods of geometrie quantization [60,61], we foUow Heisenberg's philosophy 
based on careful analysis of spectroscopic data. From his point of view we 

^**In view of Eq(60) on Page 128 of the hook by Dirac [22], the corresponding path integrals 
can now be easily constructed. In this paper for the sake of space we are not going to take 
advantage of this observation.We refer our readers to Ref. [59] for an illustrative example. 
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have the set of classical observables {Oi{t)} which is supposedly complete. This 
means that treating the Poisson brackets as Lic brackets we have 

{0.,0,} = ^C^.Ofc. (6.4) 

k 

Accordingly, quantum mcchanically, instcad of Eq.(6.3) we need to consider the 
result 

[Ò,i?]=^CÌÓfc (6.5) 
fe 

vaUd for any ti So that under such circumstances (quantum) dynamics for- 
mally disappears! This observation can be strengthened due to the foUowing 
chain of arguments. In mathematics (see Part III, Section 3.2) expression hke 
U{t)OiU~^{t) = Oi{t) = AdjjOi defines an orhit for the operator Oi in the Lie 
algebra (made of operators {0^}) caused by the action of elements C/ from the 
associated with it Lie group. At the same time, the mathematics of Lie groups 
and Lie algebras produces for [0,H]= adfjO where both O and H are in the 
Lie algebra {Oij.Evidently, we can obtain the same (or even greater) Informa- 
tion working with Ad operators instead of ad. In particular, we would like to 
consider the trace, i.e. tì:{AduOi} = x{Oi), which is just the character for Ój. 
Clearly, it is time-independent. If this is so, then, what is the meaning of an 
orbit ? This topic was discussed at length in Parts II and III of our work. To 
avoid repetitions we refer our readers to these papers. If there is no time evolu- 
tion for the character, superficially, nothing happens. This is not true, however 
as was recognized already by Dirac, Ref.[22]. In Chapter 9 he writes: " The 
Hamiltonian is a symmetrical function of the dynamic variables and thus com- 
mutes with every permutation. It foUows that each permutation is a Constant 
of motion. This happens even if the Hamiltionian is not constanil^." Hence, 
the orbit AdjyOi is caused by permutations. These can be analyzed with help 
of the torus action thus leading to the Weyl-Coxeter reflection group W — N/T 
described in Section 3.1 of Part III and to the associated with them Lie algebras 
discussed in section 3.2. of Part III. Representations of these Lie algebras (in- 
cluding the afhne Lie algebras) produce ali known quantum mechanical results 
as well as those of conformai field theories (CFT). This was explained in Part 
Illlf^ At this point we would like to do more. 

^^That is time-dependent. 

•^"Many quantum mechanical problems do involve time evolution, e.g. decay of the 
metastable state, etc. To account for such phenomena we should consider random walks 
on groups. An excellent introduction to this topic can be found in the monograph by Diaco- 
nis[62]. 
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6.3 Symmetric group and its relatives 



As is well known, the symmetric group Sn has the foUowing presentation in 
terms of generators Si and (Coxeter) relation^ 

sì = 1 
SiSj — SjSi for |i — j| > 2, 
SiSi+is, ^ Si+is,s,+i. (6.6) 

If there is a set of n elements (say, the Weyl roots arranged in a certain order) the 
generator Si interchanges an element i with i + 1 so that Si, s„_i generate 
5'n.Clearly, there are ni permutations in the set of n elements. If we assign 
the initial ordered state, then any other state can be reached by successful 
application of permutational generators to this state so that the word w = 
SaiSa2 ■ ■ ■ Sai (whcre the indices ai, a; represent a subset of the set of n — 1 
elements) can be identified with such a state. Since one can reach this state in 
many ways, it makes sense to introduce the reduced word w whose length l^w) 
is minimal. With these definitions, we would like to complicate matters a bit. 
We would like the generators of Sn to act on monomials x"* = a^^^Xj^ • • • 
FoUowing Lascoux and Sciitzenberger (L-S) , Ref. [63] , we introduce an operator 
di via rule 

a, il^ilL. (6.7) 

It acts on monomials such as x'* in such a way that the generator Si acting on 
the combination converts it info a;^'^^^^^]^. By design an action of this 

operator on monomiai is zero if = a^+i, otherwise it diminishes the degree of 
the monomiai by 1. In addition, these authors introduce the operators ttì 

TTi = Xi+1 (6.8a) 

and 



(6.8b) 



Evidently, in view of Eq.(6.8b), it is sufhcient to use just one of these operators. 
Because of this, foUowing Ref. [63], we introduce an operator Di{p, q.r) = pdi + 
qTTi + rsi with p, q, r being some numbers. L-S demonstrated that such an 
operator obeys the braid-type relations (the 2nd and third of Eq.s(6.6)) while 
the relation = 1 in Eq.(6.6) is replaced by 

DÌ qD, + r{q + r). (6.9a) 

As is well known, the last relationship (with constants q and r properly chosen) 
defines the Hecke algebra 7J„ of the symmetric group Sn- For the future use we 
shall rewrite it in the commonly used form as 

Dì = (1 - Q)A + Q- (6.9b) 



A quick introduction can be found in Appendix A of Part II. 
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Hn should be considered as a deformation of Sn- To be precise, such defined 
Hecke algebra is of the A„_i type in the Coxeter -Dynkin classification schemj^. 
Since connection of Hecke algebra with knot theory is well known [64], we 
are not discussing it in this work. Instead, we would like to connect these 
results with traditional quantum mechanics thus bringing back the spirit of old 
ideas of Maxwell, Kelvin and Tait [65]. From such point of view the differences 
between quantum mechanics, quantum field theory and string theory practically 
disappear. 

FoUowing Kirillov [58] , we begin with rclabcling prcviously defined operator 
di as bij. Next, let di be the usuai operator of differentiation, i.c. di = thcn 
we define the Dunkl operator Vi by 

V,=d,+kY,K. (6.10) 

where k is some (known) Constant. Such an operator acts on monomials (poly- 
nomials). It possess the property wViW^^ — V^i^i) \/w G Sn.Consider now the 
commutator [DijVj]. It can be rather easily demonstrated [58] that such a 
commutator is zero if bij satisfy the classical Yang-Baxter equations (CYBE) 

[bi2, bis] + [òi2, 623] + [^13, ^23] = 0. (6.11) 

Conversely, Eq.(6.11) can be taken as a definition of bij. In such a case we no 
longer need its explicit form given by Eq. (6. 7). This is facilitatcd by the designing 
of the so called degenerate affine Hecke algebra. Such an algebra is made as a 
semidirect product of 5„ with the familiar commutator algebra 

Xi+iSi - .SiXi = h, XiSj = SjXi Vi 7^ j,j + 1 (6.12) 

where h is some Constant analogous to It should be clear at this point that 
Eq.s (6.12) are discrete analogs of the Heisenberg commutaton rule, Eq.(5.8). 
Let US introduce yet another operator s; = Si + /lòj.i+i. It is designed in such a 
way that it obeys the braid relations: 

sihsi = S2S1S2. (6.13) 

Moreover, if we define R12 = siSi,iÌ23 = S2S2,^i3 = •sii?23Si = S2Ì?i2S2, then 
the above Eq.(6.13) becomes equivalent to the standard Yang-Baxter (Y-B) 
equation for Rij = 1 + hbij ( or Rij ~ exp(/iòy ) for /i — s- 0), i.e. 

R12R13R23 — -R23-Rl3-Rl2- (6.14) 

Based on this logie, it foUows, that the quantum Y-B equation, Eq.(6.14) for 
Rij implies the classical Y-B Eq.(6.11). 

•^^Because of this, one should be aware of existence of Hecke algebras for other type of 
refiection groups [66]. We are not going to use them in this work. 

'^''In fact, it is equal to h in most of cases known in hterature. In this work we do not impose 
such a requirement. 
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Ali this discussion looks a bit formai. Indccd, why to introduce the operator 
Why to be concerned about the commutator [DijVj]! What the Yang- 
Baxter equations have to do with ali results of earlier sections? We would like 
to provide answcrs to thcsc questions now and in the ncxt scction. 

First, consider an equation T>if = 0. It can be written alternatively as 



which is just the celebrated Knizhnik-Zamolodchikov (K-Z) equation. This 
means that: a) the operator Vi is effectively a covariant derivative (the Gaus- 
Manin connection [53] in the formalism of fibcr bundlcs) and b) that the van- 
ishing of the commutator [DìjVj] is just the zero curvature condition [53,67]. 
The question stili remains : how Qij in Eq.(6.15) is related to 6^ ? The answer 
was found by Belavin and Drinfeld [68] and summarized in Ref.[69], page 46. 
In the simplest "rational" case we have bij{z) = as expected. More com- 
plicated trigonometrie and elliptic cases were found in Ref.[68] and summarized 
in Ref.[69]. From these refercnccs it should be clcar that sincc solutions to 
the K-Z equations are expressible in terms of hypergeometric functions of single 
and multiple arguments, ali examples of exactly solvable quantum mechanical 
problems (including those involving the Dirac equation) found in textbooks on 
quantum mechanics are covered by the formalism just described. 

At this point it is legitimate to ask: ali this is interesting but not new. How 
these results are related to the honeycombs and puzzles discussed in earlier 
sections? We provide an answer to this question below. 

7 Back to fusion 
7.1 Motivation 

In this subscction wc would like to study the foUowing. 

Problem 7.1. To what extent the fusion rule, Eq.(3.7), valid for characters 
sx of symmetric group Sn should be modified if instead of this group we consider 
its deformation caused by our use of the Hecke algebra _ff„? 

We provide an answer to this problem having the foUowing goal in mind. 
Almost simultaneously with publications of KT honcycomb papers there ap- 
peared a publication by Gleizer and Postnikov (G-P), Rcf. [70], where graphical 
methods alternative to those developed by KT were used, essentially for the 
same purpose of calculating the L-R coefììcients. These alternative graphical 
methods involve braids, the Y-B and the tetrahedron equations. Since these 
equations play only an auxiliary role in G-P's work, many things where left 
unexplained. For instance, as soon as one introduces these equations one leaves 
the domain of symmetric group Sn and enters the domain of Hecke algebra for 
this group. From G-P work it foUows that the fusion rule, Eq.(3.7), is expected 
to remain the same. This happens to be the case most of the time but not 
always! The proof can be found in the paper by Wenzl, Ref.[71], Theorem 2.2. 




(6.15) 
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In the case if Q in Eq.(6.9) is the m-th root of unity the fusion rule, Eq.(3.7), 
should be replaced by more elaborate fusion rule to be discusseci in the subsec- 
tion 7.4. The diagrammatical methods developed in G-P work provide no clues 
regarding the possibility of such an alternative. It should be noted though that 
the KT graphical methods also fall under the same circumstances. In the fol- 
lowing subsections we provide evidence that the "anomalous" case corresponds 
to the situation when already familiar L-R coefficients should be replaced by the 
Gromov-Witten (G-W) coefficients (invariants). In some cases (to be specified) 
the fusion algebra, Eq.(3.7), is replaced by the Verlinde-type algebra. Nev- 
ertheless, since the L-R coefficients obtained with help of KT diagrammatic 
methods can be used as an input info more complicated expressions for the 
G-W invariants, e.g. read Appendix B, this justifies their place in this work. 

Remark 7.2. In view of Eq.(6.9) and the fact that Rij = f -I- ft-ò^ , it should 
be clear that representations for both the Hecke algebra and Yang-Baxter equa- 
tions are interrelated (and even coincide!). This is indeed the case as demon- 
strated by Jimbo, Ref. [72]. Alternative derivations can be found in the peda- 
gogically written paper by Ram, Ref.[73jf3. For non exceptional (generic) Q's 
calculation of characters of Hecke algebra is nicely explained in the paper by 
King and Wybourne, Ref. [75]. Since these are deformations of Schur functions 
s\{'x.) that are smoothly dependent on Q, the fusion rule, Eq.(3.7), remains 
unchanged. 

The Information just described is sufficient to bring us to our ncxt topic. 

8 Mapping class group 

To understand better what foUows, some facts about the mapping class group 
are helpful at this time. We discuss them bere using pedagogically written 
paper by Jones [76]. Consider some Riemann surface TZg of genus g. Every 
orientation-preserving homeomorphism of TZ is isotopie to the product of Dehn 
twists [77]. As is well known, e.g. see [78], every TZg admits pants decomposition 
into coUection of the trice punctured (holed) spheres. This decomposition can 
be made along Ci, c^g~i simple (non intersecting) closed curves. Every Dehn 
twist can be represented as a combination of Dehn twists around just described 
set of "basis" curves as demonstrated by Dehn. Subsequently it was realized 
that it is sufRcient to bave just 2g-|-l basic curves for this purpose [77]. The 
mapping class group Mn{g) is generated by the coUection of Dehn twists around 
these basis curves modulo twists isotopie to identity. To understand properties 
of this group it is convenient to consider TZg as branched covering (2-to-l) of 
the sphere S'^ with branching done at 2g + 2 points as depicted in Fig.l3. 

Thus, S'^ — TZg/i, where i is involution depicted in Fig.l3. Let Qg = 
{qi, ...q2g+2} be the branching set of points on S'^ whilc Qg = {(ji, ...529+2} 
the corresponding set of points on TZg. From Birman, Ref. [77], page 182, one 
can find how the Dehn twists on TZg are related to the set Q. 

34 See also [74]. 
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Figure 13: Riemann surface of genus g a) is a two-fold branched covering of a 
two sphere b) 

In particular, to each Dehn twist on TZg one associates points g2j-i and 
q2i on TZg (1 < ? < g) through which the Q-th basis curve is passing. The 
Dehn twist on TZg is projected into S*^ in the forni of the homeomorphism ojì 
(1 < i < 2g + 1) of S"^ resulting in exchange between the points q2i-i and g2i 
which leaves the rcst of points fixed. If 6i denotcs the isotopy class of the Dehn 
twist about Ci, it can be demonstrated that such class obeys the braid group 
Bn relations given by 9i9j = 9j0i if \i — j\ > 2 and 9ì9ì+\6ì = 9ì+\9ì9ì+\ 
othcrwise. In vicw of this, the mapping class group of the 2(7 + 2 punctured 
sphere 52 {2g + 2) is gencratcd by a homoniorphic image of thcsc gencrators 
which can be represented by 2(7 + 2 strings (braids) gencrating the braid group. 
A presentation for M 52 (2g + 2) is given by the Theorem 4.5. of Birman's book, 
Ref.[77]. ExpUcitly, it is given by 

uJiUjj = ujjUJi if \i — j\ > 2 , 

(a;i,...,a;2,+i)2«+2 = l, ^^"'^ 

LOl ■ ■ ■ W2gWÌg+i ■ ■ ■ LVl = 1. 

The homomorphism just describcd sends 9i to luì. In view of the involution 
depicted in Fig.l3, the kernel of this homomorphism is of order 2. We are 
interested in finding out wether the Hecke algebra -ff2g+i can be associated 
with the presentation given by Eq.(7.1). 

If Si is the generator for the symmetric group S2g+i (associated with ex- 
change of 2 points on S^) we are interested in mappings of generators Tj of 
Hecke algebra -f?2g+i into s, and s, into Wj. This happens to be possible but 
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nontrivial as discusseci in the Jones paper, Ref.[76]. Tiie nontriviality comes 
from the fact that representation T/s depends on Q so that one has to satisfy 
the constraint (wi, ...,u}2g+iY^^'^ = 1 for arbitrary Q's. Let t:y{sì) = ujì be the 
desired mapping (already from Ti to Ui) with Y indicating the Young tableaux 
associated with representation of the symmetric group. Lemma 9.2. by Jones 
allows US to make it in such a way that (irrespective to the actual vaine of Q) 
one obtains 

^^(si,...,S2<,+i)'^+' = 1 (7.2a) 

where the prime for Tr'y indicates the needed "adjustment" to make presentation 
Q-independent. Jones argues that the relation wi • • • uj2g^'2g+i • • • = 1 in 
Eq.(7.1) is equivalent to (w2, t^2g+i)^^^^ = 1- Accordingly, one obtains, 

^y(s2,...,S2<,+i)'^+' = l. (7.2b) 

Based on this observation by Jones, one can continue this downsizing process 
thus obtaining the flag of Hecke algebras Hi{Q) C H2{Q) C • • •i?2g+i(Q). The 
important theorem by Jones (to be used below) can be stated now as follows 

Theorem 7.3. (Jones [76], page 361) Let Y he a Young diagram and let 
tt'y be the corresponding representation of -629+1 designed in such a way that 
Eq.(7.2a) holds (for any Q's). Then n'y defines a representation of Mg2{2g + 2) 
via LOi — > t:'y{sì) if and only if Y is rectangular. 

In other words, one begins with the rectangular Young tableaux of size mxn 
as discussed in Section 3 and obtains ali Young tableaux which can fit into 
this rectangle by sequentially deleting boxes (one at the time).To finish this 
subsection we need to discuss difference between the mapping class group of the 
sphere S"^ and that of the disc D^. Since the disc can be viewed as a sphere 
with just one point deleted the generators luì described before can be used in 
the present case as well. Thus we obtain the foUowing theorem 

Theorem 7.4. (Birman [77], page 32). Let M be a group of automorphisms 
of 7ri(_D^ — Qn) which are induced by the homeomorphisms of — Q„ which 
keep the houndary of fixed poinwise. Then M is precisely the group 

Remark 7.5. The fundamental group 7ri(-D^ — Qn) is a free group of rank 
n. This group is made of loops anchored at some point in and such that 
each loop encloses only one puncture. In the next subsection this group will pay 
a role of monodromy group for the K-Z equation. 

Remark 7.6. In the case of a sphere there is no boundary which is fixed 
pointwise. As result, one gets relations 3 and 4 in Eq.(7.1) which are simple 

^^That is Bn is the braid group and the set Q„ is analogous to that determined immediately 
after Fig.13. 
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consequences of rotational invariance of S"^ (e.g. see Fig.13 of Ref.[77], page 
157). 

Being armed with these facts we are ready to discuss the next topic. 

8.1 Monodromy group of the Knizhnik-Zamolodchikov 
equations and the Riemann-Hilbert problem 

The above discussion about the mapping class group seems to be detached from 
the rest of this paper. We would Hke to correct this deficiency now For this 
purposc, wc nccd to rcmind to oiir rcadcrs about some basic facts from the 
theory of K-Z- type equations. According to Refs. [67,79] the K-Z Eq.(6.15) is 
equivalent to the following system of equations 

^ f) ^ fi 

^ — / = and ^iY^z, — )f = (^1^,,)/. (7.3) 

i=l ' 1=1 ' i<j 

The first of these equations comes from the fact that flij = flji. The second 
is obtaincd from Eq.(6.15) by multiplying both sidcs of this cquation by Zi, 
summing over i's and again taking into account that = iìji. The first of these 
equations indicates that solutions should be translationally invariant so that 

only diffcrcnccs of argumcnts must bc uscd. The second impHcs that solutions 
must bc homogcnous. To cxtract this homogencity we foUow Rcf. [69],page 37, 

and introduce new variables (,i = — i = 1, ...,n — 1, = Zn- In terms of 

Zi 

these variables the K-Z equations acquire the following standard form 

d 

'ii-^f = ^«("^i) •••> W-i)/ , « = 1, n, (7.4) 

whcrc the A[s are holomorphic functions in the domain D :={|<>jl < 1, 1 < J < 
n}. Such form of the K-Z equations allows us to use general locai theory of 
linear differential equations of Fuchsian-type. In particular, the theorem which 
is formulated and proved on page 121 of Ref.[80] provides the desired solution. 
It is given by 

/ = ^o(^i,...,^^)^f^°^---^^(°), (7.5) 

where Fo(<ji, <;„) is an n x n matrix- valued function holomorphic in the 
domain D := < 1,1 < j < m} and such that Fo(0) = /d.These results 
allow US to mention about the Riemann-Hilbert (R-H) problem and by doing so 
to make a connection with previous subsection. The results which foUow will 
bring US directly to the discussion of the G-W invariants. 

To discuss the R-H problem in the domain D we would like to rewrite 
Eq.s(7.4) in the following equivalent Fuchsian-type form [80] 

|-/ = ij, (7.6) 

where Ai = Aj(<ri, Cn-i)/<?i. Consider now the limiting case: Q ^ for a 
subdomain Di where presence of other singularities can be neglected. Without 
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loss of gcncrality wc can subdividc D into such subdomains so that the total 
(global) solution in D is made of locai solutions in respective subdomains. 
The requirement that these solutions must agree in the overlapping regions 
is the csscncc of the R-H problem stated in simplc tcrms. The n = 2g + 2 
punctured sphere discussed in the previous subsection can be looked upon as 
2g + 3 punctured disc. Hence, we can initially develop our discussion for the 
disc D. In both cascs: D or 5^, the motion of subdomains Di on these surfaces 
is controlied by the action of the respective braid groups as explained in the 
previous subsection. 

Consider now a monodromy (holonomy) around given singularity. The ex- 
istence of this monodromy is assured by the fact that Eq.(7.6) is linear equa- 
tion whose solution is defined with accuracy up to a Constant matrix which 
wc shall cali Ti. This matrix can be found by notieing that the matrix Ai 
has only the first order (Fuchsian) pole as singularity. Going around this pole 
once will pie up a phase resulting in the monodromy matrix Tj. It is given by 
Ti = cxp{2TTÌAi{0)). This can be easily understood with help of Eq.(7.5) taking 
into account the analytical properties of the logarithmic function. Consider now 
Eq.(7.5) in the overlap of two domains Di and Dj. Evidently, in such domains 
we can use just one locai coordinate so that at the overlap of these domains one 
has 

Fo{<;i)F,-\,j)=<;t^^%^^^°\ (7.7) 

With help of Eq.(7.7) the R-H problem can be formulated now as a problem of 
finding of a holomorphic vector function f((r) with good behavior at oo and 
such that in the complex piane it obeys an equation f_(<r) = M(^) f+(<;) for 
some prcscribed nxn matrix M{<;) and for a contour (closed) C in <^— piane such 
that f+ (?) and f_ (?) He respectively inside and outside of the domain enclosed 
byC. 

This problem can be reformulated a bit clifFerently as follows. Given a set of 
n points {ai, a„} in C<r and nxn matrices Ai, An representing the mon- 
odromy group G of these points fìnd ali equations of the type given by Eq.(7.6) 
which have the monodromy group G- We can complicate matters further by mak- 
ing the set of points to move in C? (or S^). This leads to the isomonodromic 
deformation problem. It can then be formulated as follows. 

Problem 7.7. (Isomonodromy problem) For a given representation of the 
monodromy group G find dependence of matrices A^ in Eq.(7.6) on location of 
poles given by the (moving) set {«i, ...,a„}. 

This problem was solved by Schlesinger [81]. We would like to provide some 
needed details within the context of K-Z equations. In particular, taking into 
account that Qij — flji Eq.(6.15) can be rewritten as follows 

df = Tf, where 

r = Yl '^-^{dzi-dzj). (7.8) 

l<i<j<n ■' 

The previously imposed requirement \Pi,Vj] = 0, e.g. see Eq.s(6.10),(6.11) 
can be now rewritten as d{df — Tf) = implying the Frobenius- type equation 



40 



r A r = which holds only if 



-Oforz^j^fc^Z (7.9a) 

and 

, n,k + fìjfc] = [^tj + , i^jk] ■ (7.9b) 

According to Kohno, Ref.[64], these are the infinitesimal pure braid relation^. 
For n = 3 they coincide with earlier obtained CYBE, Eq.s(6.11). Kohno, 
Ref.[64], demonstrated that, at least in the case of rational solutions of the 
CYBE's, the results, Eq.s (7.9 a,b), can be brought into correspondence with 
the CYBE's for n > 3, i.e. for 1 < ? < n. This result provides an independent 
support of earher obtained Eq.(6.11) and conneets us with the mapping class 
group presentation, Eq.(7.1). Evidently, the braid relations in Eq.(7.1) (adopted 
for iS*^) become the YBE's. Since the obtained isomorphism involves earUer dis- 
cussed monodromy matrices we come to the conclusion that the monodromy 
representation for the K-Z equation is equivalent to the YBE representation. 
This is known as the Kohno-Drinfeld theorem, Ref.[67], Thm 19.4.1. The in- 
troduced concepts even though being useful, play only an auxiliary role in this 
work. They were introduced mainly for the sake of the discussion presented in 
the next subsection. 

8.2 The multiplicative Horn problem 
8.2.1 Emergence of Gromov-Witten invariants 

In view of results we just obtained and, taking into account Eq.s(7.1), the mon- 
odromy matrices for the punctured sphere should be subjected to the fol- 
lowing constraint 

n 

J|exp(z2^A,) = 1, (7.10) 

i=l 

where I is the unit matrix and n = 2g+l. Taking into account Remark 7.5., this 
equation has a simple geometrical meaning. It represents loops (holonomies or 
monodromies) around n + 1 points, that is it represents the fundamcntal group 
TT of S*^ with points {ai, a„_|-i} deleted. As it is written, this equation suffers 
from the fact that it is not reflecting the differences between the topology of the 
disc D and that for the sphere S^. Because of the Remark 7.5., this equation 
cannot be used as such for the disc since the fundamental group Tri(D^ — Qn) is 
free group of rank n. Even though it can be used for S'^ it does not reflect the 
constraint Eq.(7.2) adequately. Fortunately, this deficiency is easily correctable 

^^The difference between the pure braid and braid groups is exactiy the same as the differ- 
ence between the statistics of distinguishable (colored) and indistiguishable (colorless) particles 
[77]. Since at the fundamental levei elementary particles, say electrons, are indistinguishable, 
in 2 dimensions their motion is described by the braid group Bn. Thus, the braid group is 
made of a semidirect product of pure braid group F„ and the pcrmutation group Sn- 
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if we rewrite Eq.(7.10) in the alternative forni as foUows 



exp(i27rAj) = exp(i27rdl) (7. Ha) 

with d — 0,1, 2, ... or, even more generally, as 

Yl exp(i27r--^) = exp(i27rl). (7.11b) 

Matrices exp{i2TrAi) are unitary by design and each of A'^s is diagonalizable 
so that X{Ai) = {Al (Ai), Xk{Ai)} represents the eigenvalue set for the matrix 
Ai. Since matrices exp(i27rAi) are unitary their determinant is 1. This leads to 
the requirement Xi{Ai) + ■ — h Xk{Ai) = 0{N moddi) Vi provided that di = 
d. For the sake of comparison with earUer sections we would Uke to consider 
(without loss of generaUty) the case of n = 3. Thcn, instead of Eq.(1.6), we 
obtain. 

Al H Afe + Mi H h Mfe = i^i H h i^fe + Ndi + Nd2 + Nds (7.12a) 

which is essentially the K-M-B-B condition, Eq.(1.5). The same resuh can be 
rewritten as 

|A| + 1^1 = + A^d. (7.12b) 

In view of this relation, the fusion rule, Eq.(3.7), should be modified accordingly. 
Following Ref . [27] , we write 

ax^a^^Y.l'C^x^.idy^- (7-13) 

d,iy 

The star symbol represents the product of "quantum" cohomology classes. For 
d = this symbol becomes again the usuai dot symbol used in Eq.(3.7). The 
"quantum" L-R coefficients C'(^{d) are in fact the gcnus zero 3-point Gromov- 
Witten invariants which are structure constants in the "small" quantum coho- 
mology ring [82]. 

The multiplication law, Eq.(7.13) takes place only if Eq.(7.12b) holds. 
This equation replaces earlier equation |A| + |/^| = used for computations 
in Eq.(3.7). The r.h.s. of Eq.(7.13) is a polynomial in q (where q stays for 
"quantum "). Its physical role is clear froni Eq.(7.13): it plays a role of fu- 
gacity associated with the degree of mappin sFI d. More generally, in view of 
Eq.(7.11b), we should replace q"^ in Eq.(7.13) by g'^ — Y\qf'- To avoid unnec- 

i 

essary complications we shall be working just with d from now on. Evidently, 
under such circumstances the q indeterminate becomes an analog of Q in the 
Hecke algebra, Eq.(6.9b), of the symmetric group. Connections with the Hecke 
algebra are highly nontrivial. They are discussed below and in Appendix C. In 



•^^Prom the punctured sphere to the Grassmannian G{m,k) as can be seen from 
Eq.s(7.2),(3.7) and (3.8) 
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order to prepare our readers for this discussion we would like to proceed with 
actual computation of the G-W invariants. 

There are many ways to compute these invariants. We would like to discuss 
only those which are logically compatible with results presented in previous 
sections. In particular, we would like to connect the results presented in Section 
5.3 with what has been discussed now. FoUowing Walton, Ref.[83], we begin 
with the Weyl character formula, our (III. 2. 28). It is given by 

X(A) = "-(^)eH (7-14) 

so that the fusion rule (Walton's Eq.(1.5)) reads: 

X(A) • x(a*) = E C'^M'')- (7-15) 

For the sake of space we refer to our Parts II (Appendix A) and Part III, Sections 
1 and 2, for ali definitions and notations. In Humphrey's book, Ref. [84], on 
page 140 our readers can find (our) Eq.(7.15) with details of its derivation. 
Some of these details are discussed below. Clearly, since the constants C^^ 
will be different for different Weyl-Coxeter groups, the formulas above include 
the fusion rule, Eq.(3.8), as special case (for characters of symmetric group 
Sn which is of the type A„_i in Dynkin's classification of the Weyl-Coxeter 
reflection groups). The Kostant multiplicity formula, Eq. (III. 2. 31), plays a very 
important role in Walton's analysis. In particular, he uses it in order to arrive 
at the Steinberg's formula for C^^ 

^Am= e £{^w)P{wX + Vfi-,y). (7.16) 

In arriving at this result Eq. (III. 2. 31) for the Kostant multiplicity formula was 
used. Detailed derivation of the Steinberg's formula can be found on page 141 
of Ref. [84] . Explicit use of this formula is inconvenient though. It is given bere 
to emphasize the combinatorial and symplectic nature of the L-R coefficients 
in accord with Section 5. 3 [fi. Since symplectic nature of C^^ was emphasized 
in Section 5.3., we brought Eq.(7.16) to the attention of our readers with 
additional purposes in mind. Specifically, Eq.(7.16) can be used not only for 
computations involving more traditional Weyl-Coxeter reflection groups but also 
for their affine extensions (e.g. see Appendix A for Part II for definitions). In 
this case Eq.(7.16) formally stays the same, except now the Weyl group W is 
replaced by the affine Weyl group VF'*^' (or M^('")) where k (or m) have the 
same meaning as in our Eq.(3.3) (or (3.4)). This fact is not sclf-obvious and 
will be explained. Ecfore doing so we note the foUowing. Using Eq. (III. 2. 16) we 
know that e(t)w) = (-l)H'"") = (-1)'(>')(-1)'("'). This observation allows us to 

■^^More details on symplectic interpretation of Eq. (7.16) can bo found either in Part III or 
in recent paper by Guillemin and Rassart, Ref. [85]. 
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Figure 14: Legai a) and illegal b) rim hooks 

rewrite the affine version of Eq.(7.16) in equivalent forni given by 

Clf= E ^HC^" (7.17) 

with C™^" representing standard L-R coefficient which can be calculated with 
help of honeycombs and puzzles as previously discussed. The fact that this is the 
case, unfortunatcly. is net sufFicicnt for the efficient calculation of C'^^\ Hcnce, 
the task lies in finding the efficient scheme for such calculations. By doing so 
connections between the coefìicients C^^^f* and C'^^{d) will become apparent. 

We begin with few definitions. In particular, in connection with tlie Young 
diagram defined by partition A, we define an n rim- hook of this partition. It is 
a connected subset of n boxes of A such that it does not contain a 2 x 2 square. 
A rim-hook is legai if by removing it from the Young diagram A the remainder 
is stiU be a vahd Young diagram. Othervise, the n rim- hook is considered to 
be illegal as depicted in Fig. 14. 

An n-core for the partition A corresponds to a partition obtained by sequen- 
tial removal of legai n-rim hooks till one obtains the configuration for which 
such removals are no longer possible. Let |A| be the weight of partition A, 
i.e. |A| = Al + A2 + • • •. Let then the weight of the resulting core partition 
be |core„A| . The number r„(A) =(|A| — |core„A|) /n is the number of n— rim 
hooks removed in this process. Furthermore, defìne the width{R-) of an i-th 
n-rim hook Ri as the number of columns it occupies in the Young diagram. 
Let w i=width{Ri) denote this number. As in Section 3.1, we place our originai 
Young diagram into m x k rectangle so that m + k = N. Because of this, we 
write C^^id) = C'(^{d-,m,k) so that if, say, A is the partition conjugate to A, 
one can show [86] that 

C^^(rf;m,fc) = Cf_(d;/c,m) (7.18) 

in accord with the fact that Cj^^ = C^^ for the ordinary L-R coefRcients. To 
make actual computations using these definitions we need to use the foUowing 
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Lemma 7.8. (Bcrtram, Ciocan-Fontaninc, Fulton [86], page 733). // ^ is 
the partition resulting from removing an n-rim hook from X, then the (quantum) 
cohomology class a\ is related to the (classical) cohomology class through 
the following relation 

<Tx = (-l)'=-'"<za^ (7.19a) 

where w is the width of the n-rim hook and k is the width of the Young tableaux 
containing partition X.If A contains an illegal n-rim hook or if A^+i > and 
A contains no n-rim hooks or if Ai > k, then ax = 0. 

The abovc lemma allows us to use previously discusseci fusion rule, Eq.(3.7), 
for calculation of the Gromov-Witten coefficients. We would like to show step-by 
-step how this is done. Let both partitions A and fx belong to the mxk rectangle 
and i?r„(A) be the respective n— rim hooks removed, then, instead of 

Eq.(7.19a), we obtain: 

ax = e{X/,j)q^-^^^a^, , (7.19b) 

r„(A) 

where e(A//i) = Yl (^1)*^ If the partition fi is net contained in m x fc 

i=l 

rectangle, then = 0. Next, we combine the fusion rule, Eq.(3.7), with just 
obtained result in order to obtain the prescription for calculation of the G-W 
invariants. Appcndix B contains an cxample of calculation of G-W invariants 
in which the KT honeycombs are used for calculations of C^^ as an input. In 
actual illustrative calculations done in Appendix B we took into account that 
d = rjv(A) and n ^ N in accord with Eq.(7.12b). Thus, in general, we obtain: 

o-A * cr^ = J2 E éCx^{d; m, k)cj^, where Cj^^(rf; m,k) = J2 ^{p/^)Cxu,- 

i^gfexmd=0 p 

(7.20) 

This result admits a somewhat different interpretation. For instance, let us 
introduce the " quantum" cohomology class (jp via 

^P = E q'' E <pl^)'y., (7.21) 

then we can formally rewrite Eq.(7.20) as 

c7x*a^ = j:C'x^cTp, (7.22) 
p 

where C^^ is the classical L-R coefHcient computable with help of KT honey- 
combs, puzzles or hives. This result is only in formai agreement with earlier 
obtained Eq.(3.7) since in Eq.(3.7) there is no restrictions on summation over p 
while in the present case C^^ = C^^{m,k) and is zero otherwise. Significance 
of such a restriction is discussed in the Appendix C. Obtained result is easy to 
understand using physical arguments. Indeed, using Eq.(7.12b) we notice that 
with N ^ oo the only way to satisfy this equation is to let d = 0. Hence, in view 
of the Thcorcm 5.1. Eq.(7.22) makes sense. Clearly, it should be understood 
only as qualitative result since summation over p is actually restricted. 

We have accumulatcd cnough rcsults cnabling us to inject more physics into 
what was obtained thus far. This is done in the next subsection. 
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8.2.2 Verlinde algebra and Hecke algebra at the root of unity 

Whilc the mathematica! mcaning of G-W invariants, cspccially for small quan- 
tum cohomology ring, is discusseci in many places, e.g. see Refs. [82,87,88], to 
our knowledge, tliis literature does not contain any information about physical 
significance of thcsc invariants. Wc would likc to corrcct tliis dcficicncy. By 
doing so we also will be able to sketch some missing links between these topics 
and those discussed in previous sections. We hope that our readers will use this 
material along with that of Appendix C as a point of departure for much deeper 
and thorough study. 

We begin with observation that there are two ways to deflne the L-R coef- 
ficicnts: onc is through the composition (fusion) law, Eq.(3.8), whilc the other 
through the Fourier series-type expansion of the skew Schur function Sxj^ [37]. 
Specifically, 

SA/^ = E^?>- (7-23) 

In view of Eq.(7.22) it is only naturai to anticipate that there must be a "quan- 
tum" analog of Eq.(7.23). Such an analog was indeed recently foimd by Post- 
nikov [27]. We would like to connect his results with those in Ref.[86] using 
some results from the paper by McNamara [89]. This will allow us to inject some 
physics into our discussion. 

Recali that if A and fi are two partitions such that /x C A if /Uj < Aj Vi, then 
a pair {fi, A) is called skew partition. For such a partition the Young diagram 
is madc of diagram for A with /i removed. Traditionally used notation for such 
obtained diagram is given by A//Lt = (Ai, Afc)/(/Xi, Postnikov replaced 

the m X k rectangle used in calculations in Appendix B by the torus obtained 
from this rectangle by Identification of its boundaries in a usuai way. He consid- 
ered skew partitions on such a torus and proved that Eq.(7.23), when adapted to 
this toroidal topology produces C'^^{d; m, k) as required. His derivation of this 
result makes the meaning of parameter d cspccially transparent and suitable 
for potential physical applications (to be discussed below). 

Remark 7.9. In view of Eq.s (7.11), (7.12) replacement of the rectangle 
by torus is completely naturai. Moreover, in view of the fact that in the most 
general case «z'' = ]~[ qf' , it is surprising that instead of using the multidimen- 

i 

sional torus the two dimensionai torus is sufficient for reproduction of the G-W 
invariants. This peculiarity is explained below and in the Appendix C. 

Following McNamara [89] , for illustrative purposes we choose the basic rect- 
angle R with parameters fc = 3 and m = 4. The /ti partition is chosen to be 
fj. = (2, 1) and the originai A partition is chosen to be A = (4, 4, 4, 4, 2, 1, 1). Us- 
ing the French style of writing of the Young tableaux (that is from the bottom- 
up) the skew tableaux X/fj, is depicted in Fig.l5a). Multiple copies of the same 
skew tableaux are depicted in the universal cover of the torus made of the basic 
rectangle in Fig.l5b). To recognize them in such setting one should pay atten- 
tion to the vertical lines denoted as V and the horizontal line denoted as H. 
The lower left corner of the originai rectangle R coincides with the intersection 
of H and V lines. After that, the partition n = (2, 1) to be removed is easily 



46 



I 

I 




Figure 15: The skew tableaux a) and its torus equivalent b) 

recognizable. The boxes labeled by x help us to identify the muhiple copies of 
R so that the partition A can be readily identified 

Next, we remove a shaded boxes froni A which is equivalent of removal of one 
7-rim hook (since in this example k + m = 7). This leaves us with the partition 
a' = (4,4,4,1). As it foUows from the Fig.15, such a partition stili does not 
fit the R so that we bave to remove yet another 7-rim hook (just like in the 
Appendix B but in reverse). The resulting (core) partition A" = (3,3) does fit 
the rectangle R. Such obtained skew partition can be encoded by (3,3)/2/ (2,1) 
or, more generally, X/d/ ji = (Ai, Xk)/d/ (/i^, /i;). Hence, d — 2. 

With such defincd toric skew partitions one can replace Eq.(7.23) by 

«A/d/M^ E C^.(d;m,fc)s,. (7.24) 

If in the d— sum in Eq.(7.20) only one term is nonzero, then C^jj{d;m, k) in 
Eq.(7.24) is the desired G-W invariant. In Ref.[90], page 379, Witten argued 
that for dimensionai reasons in the d— sum in Eq.(7.20) only one term should be 
nonzero. If this is the case, then C^^{d] m, k) in Eq.(7.24) is the G-W invariant 
and coincides with that given in Eq.(7.20). Moreover, if in Eq.(7.20) we formally 
choose q — 1 thus obtained fusion rule coincides exactly with that for the Hecke 
algebra, Eq.(6.9b), for which Q is the m-th root of 1, i.e. Q™ = l(recall that 
N — m + k) as it was demonstrated by Goodman and Wenzl [91]. These authors 
noticed that under such circumstances the fusion coefRcients coincide with those 
for the Verlinde algebra for su{k) Wess-Zumino- Witten (WZW) CFT model at 
the level tt^. 

^^In view of symmetry, Eq.(7.18), evidently, one can talk as well about sii(m)WZW model 



47 



Remark.7.10. Sincc conncctions with rcprcscntations of Hcckc algebra are 
insured by previous stated Theorem 7.3. by Jones, while connections between 
the K-Z type equations and theWZW-type models were established by Knizh- 
nik and Zamolodchikov [92], this cxplains (in vicw of Goodman- Wcnzl rcsults) 
why in calculations of G-W invariante presented above and in Appendix B the 
rectangular boxes are used. 

Remark 7.11. Thcsc facts are not snfficicnt for cxplanation of the mathe- 
matica! meaning of the rim-hook removals (additions) in such calculations. This 
deficiency is corrected to a some extent in Appendix C. 

Obtaincd results give us an opportunity to discuss some additional physical 
applications. 

8.2.3 Possible applications to solid state physics 

In Part I, Eq. (1,3.22), or in Eq.(1.5) of this work. wc noticcd that such an 
equation can be interpreted within a context of solid state physics. In view 
of earlier obtained results we would like to discuss this connection in some 
detail now. In solid state physics the Bloch thcorcm [5] rcflccts the diffcrence 
between the wavc fimction in the vacuum and that in the periodic solid. Account 
for periodicity rcquires the wave function to be written in the form \I/k(r) = 
e'^^''^Uk{r), whcre Uk(r) is some periodic function. The vector k is determined 
by the condition 

e^-^ = 1, (7.25) 

where 1 is the minimal translation vector of the direct lattice. Since this type of 
equation wc havc cncountered already in Part II, e.g. see Eq.(II.9.18), we know 
that solution is given by 

k-l = 0mod27rd (7.26) 

which implies that the vector k should belong to the reciprocai ( or dual) lattice 
and d should bc some intcgcr. In Part I wc callcd this equation as the Kac- 
Moody-Bloch-Bragg condition (K-M-B-B condition). Glearly, Eq.(7.11) is of the 
same type. This fact underscores its solid state physics relevance. In particular, 
it causcs ali physically intcrcsting quantities to changc accordingly so that, for 
instance, the momentum vectors for electrons are defined only with accuracy up 
to some vectors of the reciprocai lattice. This causes energy of the electron not 
to be well defined (as cxplained below). Eq.(7.26) suggcsts strong links with 
number theory. The book by Terras, Ref.[93], especially Ghr.lO, provides an 
excellent starting point for development of theory of electronic band structures 
of solids (and also of molcculcs) using mcthods of number theory. Glearly, in 
doing so one has to use group theory in the number fields of characteristic other 
than zero as explained in Appendix C. 

Some uscs of the WZW-type models in solid state physics are mentioned in 
the book by Tsvelik [94]. Thus far they are restricted only to 1-1-1 dimensionai 
models whose exact solutions can be obtained by other methods anyway. We 
would like to argue that aJJ many -body problems of solid state physics in which 

at the level k. 
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lattice pcriodicity cannot be disregarded should employ mathematical formalism 
discussed in this paper. 

That this is the case can be seen from the simplest example. Indeed, 
the existing band theory of sohds [5] is essentially the theory onc electron in 
the periodic 3d lattices, e.g. ah known lattices £ are made as direct sum 
C={Z /m\Z) ® (Z/m2Z) (g) [Z/m^T?) , where mi, m2, are respective lattice 
periods. AH such lattices have 2 kinds of synimetry: point-like and spatial. The 
affine Weyl-Coxeter reflection groups are capable of describing these symme- 
tries. The word "affine" accounts for effects of translational symmetry of the 
lattice (as explained in Appendix A to Part II). This symmetry is so important 
that the presence of potential in the Schròdinger equation can be often ignored 
[95] . The band theory developed for the " empty" lattice sometimes is sufficient 
for good reproduction of cxperimental data. For the empty lattice the energy 
f (k) of a single electron is known to be [96], Chrl., 

Sil.) = ^\l. + Kf, (7.27) 

where K is any vector of the reciprocai lattice. In view of Eq.(7.25) the vector 
k can have only finite number of values [96] . Unlike the traditional band theory 
of solids dealing with one electron in periodic lattice C, wc would like to place 
several electrons into such a lattice so that the total energy of noninteracting 
(except, due to the Pauli principle) is given by 

^^;^=E|ki + K|^ (7.28) 

Since the vector K is arbitrary we choose it to be the same for ali electrons. If 
this is the case, the above expression can be rewritten as 

^E^ = E(k? + 2k,.K + K^). (7.29) 

Since the total momentum of such system of electrons should be conserved this 
requires us to write 

Eki = 0. (7.30) 

i 

Finally, summing over ali energy levels (whose number is finite), that is taking 
a trace of the corresponding matrix operator, brings us back to the equation 
analogous to Eq.(7.12a). The Pauli exclusion principle requires the total wave 
function to be antisymmetric. This requirement is satisfied by the factorized 
wave function made of product of spin and coordinate-dependent part. In the 
absence of lattice periodicity the procedure of constructing such antisymmetric 
total wave function for several electrons can be found in the book by Messiah, 
Ref. [97] . This procedure uses essentially the representation theory of symmetric 
group Sn- Presence of lattice periodicity replaces Sn by its affine analog as 
discussed in Appendix C. The above picture can be complicated by accounting 
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for cffccts of the Constant magnctic ficld acting on clcctrons in pcriodic lat- 
tices. Rigorous mathematica! treatment of this problem has been initiated in 
Works by Novikov and his coUaborators [98]. The small quantum cohomology 
ring discussed in earher subsections becomes the Novikov ring under present 
circumstances [82]. 

Acknowledgement The author would Ukc to thank Alien Knutson (U of 
California, Berkeley) for his kind permission to reproduce some figures from his 
papers (with T.Tao and C.Woodward) in this paper. 



Appendix A. Details of Heisenberg's derivatìon of the commutator 
identity [x,p\ = ih 

In this Appendix wc would likc to previde some details of Heisenberg's rea- 
soning leading to the discovery of [x,p] = ih. We do this for several reasons. 
First, this would be unnecessary should his originai Nobel Prize winning paper, 
Rcf. [18], contain ali details and would bc frce of typographical errors. Second, 
although in Section 5 we mentioned Dirac's acknowledgement (on page 177 of 
his book, Ref.[22]) of the fact that light scattering experiments associated with 
mcasurcmcnt of the rcfractivc index (or diclcctric Constant) and thcir thcorcti- 
cal interpretation had lead Heisenberg to his discovery of quantum mechanics, 
nowhere in the existing literature on quantum mechanics were we able to find 
cxposition using this historic fact as the starting point for the development of 
quantum mechanical formalism. 

At the classical level consider a gas of noninteracting atoms, better just one 
atom containing N electrons which are assumed to scattcr light independently. 
The interaction between the incoming light and such an electron is described 
with help of the combination d= /3E where d is the dipole moment of the 
electron in the atom, E is the strength of the cxternal electric field which is 
assumed to be time-dependent, and /? is the polarization tensor (in the simplest 
case it is assumed to be a scalar). In the medium the strength of the electric 
field changes as compared to the vacuum. By denoting it as D it is known that 
D=E +47rP where P=Nd. Since, at the same time, by definition, d=er we 
have to have an equation for r. It is given by 

i= + a;2r + 7r =^E(t), (A.l) 

whcrc e is clcctron's chargc and m is its mass. In writing this equation it is 
assumed that our electron is bound harmonically (with the basic frequency coq) 
and that the friction is of known (electromagnetic) nature and is assumed to 
be small. Using Fourier decomposition of r(t) we obtain, 

ri^) = -- • (A.2) 
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This equation allows us to obtain P and, hence, D as foUows 



1 

D = E + 47rP = (1 + 47riV ■. — )E = £(w)E. (A.3) 

This equation defines a complcx frcqucncy-dcpcndcnt diclcctric Constant £{ijj). 
From electrodynamics it can be equivalently rewritten as e(uj) = {n(uj) — ÌK{uj)Y 
where n(a;) is the refractive index while h{ijo) is the coefììcient of absorption. 
Using these facts we can write approximately 

1 

n(a;) = l + 2^iV ^ -, — . (A.4) 

By ignoring friction in the high frequency hmit we obtain, 

n(w) = l-27rA^ -. (A.5.) 

To account for quantum mechanical efFects, Thomas, Reich and Kuhn in 1925 
(just before the quantum mechanics was born !) have suggested to replace 
Eq.(A.4) by 

n{uj) = l + 2nN-y , (A.6) 

where, foUowing these authors, we ignored friction and introduced the oscillator 
strength fi.To reconcile Eq.(A.6) with (A. 5) we have to require X^/j = 1- This 

i 

requirement is known as the sum rule. These facts were known to Kramers and 
Heisenberg, Rcf.[99], whcrc our rcadcrs can find additional dctails. To makc our 
point and to save space, we would like to reobtain the result, Eq.(A.6), quan- 
tum mechanically using modem formalism. We refer our reader to Ref.[100], 
pages 316-319], for additional dctails. Basically, we need to calculate quantum 
mechanically the dipole moment d, that is 

i;*^evi,^dh. (A.7.) 

In this expression the wave function ìp^ is calculated with help of the stationary 
perturbation theory with accuracy up to the first order in perturbation (which 
is er-E). A short calculation produces the following result for the oscillator 
strength: 

f,^ = ^I!}^\{k\x\m)\\ (A.8) 
This result can be equivalently rewritten as 

fkm = "^^^^"' {(fc I X I m)* {k \ x\m) + {k\x \ m)* {k\x \ m)}. (A.9) 

Since, however, 

imujkm {k \ X \ m) = {k \ px \ m) (A. 10) 
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Figure 16: pi with its rim-hooks and phase factors a), the same for p2 b) 

we can rewrite Eq.(A.9) as 

fkm = \ X \ k) {k\p:;o \ m) - {m\p:;o \ k) {k \ x \ m)} (A.ll) 

in 

sirice Uhm = —i^mk- Finally, we have to require J2fkm = 1- This is possible 

k 

only if 

— (m I - PxX I m) = 1, (A.12) 
in 

QED. 

Appendix B. An example of detailed computation of C'^'i^{m, k). 

In this appendix wc would like to work out an example of computation 
of the Gromov-Witten invariant C'^'^{m, k) based on Example 1 taken from 
Ref. [86] . Our calculations differ however from those in Ref. [86] since we use the 
KT schemc for computation of the classical L-R coefficients. 

In Example 1 the basic Young tableaux rectangle is taken as m x fc = 5 x 5. 
From here we obtain: m + k = N = 10 Hencc, the Icngth of the rim-hooks 
to be used is 10. Next, we are given partitions A = (5, 4, 4, 2, 2), /i = (3,2. 1) 
and v = (2,1). Based on these data, we can calculate the weights. Thcsc are 
respectively |A| = 17, = 6 and \ì^\ = 3. Since we know already that N = 10, 
the fundamental relation, Eq.(7.11b), that is |A| + = + dN now leaves 
US no choice for d. We obtain, d — 2. But this number gives us the number 
of the 10-rim hooks to be used in our computations, that is 2. If we choose v 
as the core partition then, indeed, only 2 rim hooks will fili in the rectangle as 
depicted in Fig.l6 below 

The filliiig of the rectangle (in the present case- the square) is made in such 
a way that the added rim hooks must not be added beyond the pre assigned 
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width of the rcctanglc but, at the same time, they are allowed to occupy more 
height space than the rectangle can provide. This results in two partitions pi = 
(5, 5, 4, 3, 2, 2, 2) and = (5, 4, 4, 3, 2, 2, 2, 1) which can be read off directly 
from Fig.l6. The summation over p in Eq.(7.20) in the present case takcs place 
over Pi and P2- The phase factor e{p/v) under this summation can be easily 
computed based on the information given in Fig.l6. For this purpose we need 
to calculate the widths of the 10-rim hooks in both tableaux depictcd in Fig.l6. 
For the first tableaux the widths are wi = ò and ■W2 = ^ respectively. Hence, 
the overall phase factor e{pi/u) = (—1)^"^ • (—1)^"^ = 1. Analogously, for the 
second tableaux the phase factor e(/32/;^) = (— 1)'"'^''' • (— 1)^"'' = —l.To complete 
our calculation, we need the actual values for the classical L-R coefìicients Cjil 
and Cj^l. We calculate them graphically with help of the KT method using Fig.6 
as an examplc. This example indicatcs that we have to make some adjustments 
in the initial data for partitions in order to be able to use the KT scheme. 

Hence, let us begin with calculations of C^l. The enlarged partitions hav- 
ing the same number of cntries are A = (5, 4, 4, 2. 2, 0, 0), /i = (3,2,1,0,0,0,0) 
and Pi = (5,5,4,3,2,2,2). It should be noted that the locations of the added 
zeros are quite arbitrary and that their redundancy should be discarded in ac- 
tual calculations (that is this redundancy should not affect the numerical value 
of C^l). With these remarks we need to make a set of Y-shaped tripods of 
the type depicted in Fig.l. By doing so we need to match together the num- 
bers from A and fi partitions in such a way that the partition pi is obtained. 
Also, we must pay attention to the order in which the labeling of the honey- 
comb is made, e.g. see Fig.3. After that, the number of possibilities for such 
type of matching is Cj^l. If we disrcgard rcdimdancy of added zeros, thcn we 
straightforwardly obtain C^l = 2 . To repeat this procedure for p2 we need to 
use the foUowing partitions A = (5, 4, 4, 2, 2, 0, 0, 0), p, = (3, 2, 1, 0, 0, 0, 0, 0) and 
P2 = (5, 4, 4, 3, 2, 2, 2, 1). The result of matching now produces Cjll = 1. Hence, 
altogether we obtain: C^^{d] m, fc) = 2 — 1 = 1, in accord with Ref.[86]. 

Appendix C. Hecke algebra and Kashiwara crystals 

The purpose of this Appendix is only to provide a list of kcy rcferences 
needed in support of results of the main text. Obviously, the choice of references 
is subjective. Nevertheless, it is hoped that it provides a suffìcient background 
level needed for reading any other literature on these subjects. 

We begin with Eq.(7.22). This result can be looked upon using theory of k- 
restricted (bounded) partitions recently developed by Lapointe and Morse [101]. 
These authors were looking at the foUowing problemi how fusion, Eq.(3.8), 
is going to change if instead of standard Schur functions sx one would use 
fc— bounded Schur functions s^^^ indexed by partitions A whose first part (i.e. Ai) 
is no larger than a fixed integer fc > 1? For such functions the L-R coefìicient C^^ 

should be replaced by C'^'j^ . These authors managed to prove (see their Theorem 
18) that, in accord with our Eq.(7.22), for the appropriately chosen k (actually 
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ior k = N — 1) the cocfScicnt C^'i^^^ = C'^'^{m, k) (with N = k + m). These 
results can be looked upon from broader perspective [102]. To this purpose we 
introduce the affine symmetric group Sk generated by k elements sq, si, , .., Sk-i 
and satisfying the affine Coxeter relations 

= 1 

SiSj = SjSi , i — j = ilmodA;, (C-I) 

The generators ài should be understood as Sjmodfe if « > A;. The usuai symmetric 
group Sk embeds into Sk as a subgroup generated by so,si, , ..,Sk-i but this 
embedding is not unique: thcrc are many embeddings among which one has to 
choose the most convenient [102]. The question arises: how s^^^ is connected 
with representations for 5'^? We would like to give an answer in several steps. 
First, we recali that the Hecke algebra iJ„ is the deformation of Sn with the de- 
formation parameter Q. Characters of symmetric group sx and those for Hecke 
algebra have the same fusion algebra controlied by the L-R coeffìcients (as dis- 
cussed in the main text), except for the case when Q is a nontrivial root of 1. 
In the fast case the fusion algebra for the Hecke characters coincides with the 
Verlinde algebra as demonstrated by Goodman and Wenzl [91]. Lapointe and 
Morse [101] argue that thè fusion algebra for fc— bounded Schur functions s^^^ 
coincides with the Verlinde algebra for the WZW models. Since generators Di 
in Eq.(6.9) for the Hecke algebra are deformations of the generators for 
symmetric group Sn, one can think about deformation of generators Sj for the 
affine symmetric group Sn- In Ref.[103] one can find that such a deformation is 
described by the affine Hecke algebra Hn{Q) which is the Weyl reflection group 
for the affine Lie algebra sl„. Finally, according to Ref.s[104], [105] represen- 
tations of sl„ and those for the Hecke algebra at root of unity are interrelated 
albeit in a very nontrivial way (as explained below). Hence, s^^^ (up to a 

Constant) may coincide with representation for the Hecke algebra at y/l. s^^^ 
form a Schubert basis for the cohomology ring of the affine Grassmannian as 
it is explained by Lam [106]. Affine Grassmannian was recently discussed in 
Works by Kapustin and Witten [107] and also, independently, by Frenkel [108]. 
In both cases it was discussed in connection with applications of methods of 
number theory to string and CFT. 

Because of this, we would like to explain how the number theory enters into 
previous discussion. For this purpose, we would like to discuss the notion of 
a crystal and a crystal base foUowing Ref. [109]. A crystal is made of a set B 
endowed with the maps èi, fi : BU{0} BU{0} {i G /), where / is an index set. 
These maps satisfy the foUowing conditions: a) ejO = /jO = 0; b) éfb = f^b = 
yb € B and i G I , e) 3b,b' G B and i G I such that b' = fib if and only if 
b = èib'. 

The above definition is too abstract to see the connection with crystals. This 

deficiency is easily correctable. Using results and notations of Appendix A (Part 
II) we introduce the weight lattice (typical for the affine reflection groups) 
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by 

n li 

Pi^ = {E aitoi I ai > 0, ai e Z Vi = 1 n; E «i = ^} (C.2) 

i=0 i=0 

SO that B = _P+ = U/>o^i^- ^^"^ sanie time, Eq.(C.2) has the mimber- 
theoretic meaning (e.g.see Ref.[110], Appendix) in the case if G F, = Z/gZ 

n 

with q being some prime number. In this case the expansion P — ^ Qìujì can be 

interpreted as an expansion of a number /? which belongs to the field extension of 
the number field F,. Since such an extension corresponds to a particular solution 

n 

of the familiar (by now) equation ^ Ui = l, different solutions of this equation 

represent different numbers /3. This observation is sufficient for explanation of 
a role of the operators èj and fi in the emerging picture. Indeed, by analogy 

n 

with Section 3 we can associate with each solution of ^ aj = l the Young 

diagram. The new element in doing so lies in the fact that boxes in such Young 
tableaux should be filled with numbers which belong to the field Fg. This is 
rather easy to do, e.g. see Ref.[105]. With this Information, a crystal graph can 
be constructed as foUows. We take a set of /3's as vertices of the graph (so that 
P's ^ b's) then, the operators fi connect vertices related by b' = fib. The index 
i represents a color so that the crystal graph is directed and colored by colors 
from the set /. The connections with rcprcscntation thcory and fusion can be 
figured out now based on the observation that different /3's are in one-to one 
correspondence with different Young diagrams. Each such diagram encodes a 
character for the respective group, e.g. in characteristic zero characters for S„ 
are the Schur functions SA,etc.Hence, as we just explained, the crystal graphs 
can be used instead of KT honeycombs (or puzzles) for calculation of fusion 
coefB-cicnts [111]. 

Finally, we would like to provide a few additional details regarding the actual 
role being played by the operators fi and èi in order to connect the results of 
this appendix to that of Appendix B and Section 7.4.2. In these sections we've 
noticed that different Young diagrams can be obtained from the set of core 
Young diagrams by adding(removing) the appropriate rim-hooks. The core 
Young diagrams can be looked upon as representing a vacuum state while the 
addition(removal) of the particular rim-hook is done by the raising (lowering) 
operator. In the simplest case we are dealing with addition(removal) of just one 
box of the Young diagram. Hence, we can associate the addition (removal) of 
just one box with the operator èi (with the operator fi). Clearly, if we need to 
remove(add) a rim-hook the appropriate combination of such raising (lowering) 
operators should be used. Interestingly enough, a systematic devclopment of the 
formalism we just sketched leads to the Heisenberg commutation rule, Eq.(5.8). 
Details can be found in Ref.[112]. The same formalism ior q — 1 is used for 
description of representations of the Hecke algebra 11^(^/1), Ref.[105], page 
279. 

Remark. We would like to mention that the concept of a crystal was devel- 
oped by Kashiwara, e.g. see Ref.[113], in connection with obtaining the exact 
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solution of XXZ model in the thermodynamic limit (that is in the limit of infi- 
nite spin chain). The same authors notice that the whole formalism of quantum 
groups, vertex operators, Virasoro algebra, etc. was developed and makes sense 
either for infinite chains or for finite chains with appropriately chosen bound- 
ary conditions. In the case of finite chains without specially chosen boundary 
conditions the whole apparatus of this Appendix fails. This observation is very 
important for development of our new string- theoretic formalism. It is totally 
consistent with results presented in Parts I-III and will be further discussed 
in Part IV, especially in connection to and comparison with recent alternative 
string-theoret ic developments [114]. 
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